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functions for sparticle masses, and show that they may be significantly lighter in the pMSSMlO than in the other 
models, e.g., the gluino may be as light as ~ 1250 GeV at the 68% CL, and squarks, stops, electroweak gauginos 
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The quest for supersymmetry (SUSY) has been 
among the principal objectives of the ATLAS 
and CMS experiments during Run 1 of the 
Large Hadron Collider (LHC). However, despite 
searches in many production and decay channels, 
no significant signals have been observed Oil- 
These negative results impose strong constraints 
on i?-conserving SUSY models, in particular, 
which are also constrained by measurements of 
the mass and other properties of the Higgs bo¬ 
son , by precision measurements of rare decays 
such as Bg —>■ and other measure¬ 

ments. Overall, these constraints tend to reduce 
the capacity of SUSY models to alleviate the hier¬ 
archy problem. However, their impact on a pos¬ 
sible resolution of the discrepancy between the 
experimental measurement of (g — 2)^ and theo¬ 
retical calculations in the Standard Model (SM) 
depends on further assumptions as will be dis¬ 
cussed below. 

There have been many analyses that combine 
these constraints in global statistical fits within 
specific SUSY models based on the minimal su¬ 
persymmetric extension of the Standard Model 
(MSSM) [^. Many of these analyses assume 
that the low-energy soft SUSY-breaking param¬ 
eters of the MSSM may be extrapolated using 
the renormalization-group equations (RGBs) up 
to some grand unified theory (GUT) scale, where 
they are postulated to satisfy some universality 
conditions. Examples of such models include the 
constrained MSSM (GMSSM) [9-11, in which 
the soft SUSY-breaking mass parameters uiq and 
mi /2 are assumed to be universal at the GUT 
scale, as are the trilinear parameters Aq. Other 
examples include models that relax the universal¬ 
ity assumptions for the soft SUSY-breaking con¬ 


tributions to the Higgs masses, the NUHMl 12 


and NUHM2 (see also, e.g.. Ref. 11 ), but re¬ 
tain universality for the slepton, squark and gaug- 
ino masses. Such models are particularly severely 
constrained by the LHC searches for colored spar¬ 
ticles, the squarks and gluino, which also place 
indirect limits on the masses of sleptons and elec¬ 
troweak gauginos and higgsinos via the GUT scale 
constraints, while the direct search limits on these 


2 



















3 


particles have much less impact. 

An alternative approach is to make no assump¬ 
tion about the RGE extrapolation to very high 
energies, but take a purely phenomenological ap¬ 
proach in which the soft SUSY-breaking parame¬ 
ters are specified at low energies, and are not re¬ 
quired to be universal at any input scale, a class 
of models referred to as the phenomenological 
MSSM with n free parameters (pMSSMn) [Id] . 
This is the framework explored in this paper. 
Favoured mass patterns in a pMSSMn analysis 
might then give hints for (alternative) GUT-scale 
scenarios. 

In the absence of any assumptions, the pMSSM 
has so many parameters that a thorough anal¬ 
ysis of its multi-dimensional parameter space is 
computationally prohibitive. Here we restrict 
our attention to a ten-dimensional version, the 
pMSSMlO, in which the following assumptions 
are made. Motivated by the absence of signifi¬ 
cant flavor-changing neutral interactions (FGNI) 
beyond those in the Standard Model (SM), we as¬ 
sume that the soft SUSY-breaking contributions 
to the masses of the squarks of the first two gen¬ 
erations are equal, which we also assume for the 
three generations of sleptons. The FCNI argu¬ 
ment does not motivate any relation between the 
soft SUSY-breaking contributions to the masses 
of left- and right-handed sfermions, but here we 
assume for simplicity that they are equal. As a 
result, we consider the following 10 parameters in 
our analysis (where “mass” is here used as a syn¬ 
onym for a soft SUSY-breaking parameter, and 
the gaugino masses and trilinear couplings are 
taken to be real): 

3 gaugino masses : Mi_ 2,3 , 

2 squark masses : = niq^ ^ rriq ^, 

1 slepton mass : , 

1 trilinear coupling : A, (1) 

Higgs mixing parameter : /r , 

Pseudoscalar Higgs mass : , 

Ratio of vevs : tan f3 . 

All of these parameters are specified at a low 
renormalisation scale, the mean scalar top mass 
scale, Msusy = dose to that of elec- 


troweak symmetry breaking. 

In any pMSSM scenario such as this, the dis¬ 
connect between the different gaugino masses al¬ 
lows, for example, the U(l) and SU(2) gauginos to 
be much lighter than is possible in GUT-universal 
models, where their masses are related to the 
gluino mass and hence constrained by gluino 
searches at the LHG. Likewise, the disconnect be¬ 
tween the different squark masses opens up more 
possibilities for light stops, and the disconnect be¬ 
tween squark and slepton masses largely frees the 
latter from LHG constraints. 

An important feature of our global analysis is 
that the possibilities for light electroweak gaug¬ 
inos and sleptons reopen an opportunity for an 
significant SUSY contribution to {g — 2)^ in the 
pMSSM, a possibility that is precluded in sim¬ 
ple GUT-universal models such as the CMSSM, 
NUHMl and NUHM2 by the LHC searches for 
strongly-interacting sparticles. As we discuss in 
detail in this paper, the pMSSMlO’s flexibility re¬ 
moves the tension between LHC constraints and 
the measured value of {g — 2)^ [^, with the re¬ 
sult that the best fit in the pMSSM 10 has a global 
probability that is considerably better than in 


the CMSSM, NUHMl, NUHM2 or SM. 

The main challenges for a global fit of the 
pMSSMlO are the efficient sampling of the ten¬ 
dimensional parameter space and the accurate 
implementation of the various SUSY searches by 
ATLAS and CMS. As in 16 , here we use the sam¬ 
pling algorithm MultiNest 17 to scan efficiently 
the pMSSMlO parameter space. To achieve suf¬ 
ficient coverage of the relevant parameter space, 
approximately 1.2 x 10® pMSSMlO points were 
sampled. However, confronting all these sam¬ 
ple points individually with all relevant collider 
searches is computationally impossible. In order 
to overcome this problem and still to apply the 
SUSY searches in a consistent and precise man¬ 
ner, we split the LHC searches into three cate¬ 
gories. In the first category we consider inclu¬ 
sive SUSY searches that mainly constrain the pro¬ 
duction of coloured sparticles, namely the gluino 
and squarks. To apply these searches to the 
pMSSMlO parameter space, we follow closely an 
approach proposed in 18 , which uses a vari¬ 


ety of inclusive SUSY searches covering differ- 
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ent final states to establish a simple but accurate 
look-up table that depends only on the gluino, 
squark and LSP masses. Then, in order to im¬ 
plement the other two categories of LHC con¬ 
straints on the SUSY electroweak sector and com¬ 
pressed stop spectra, we treat the LHC searches 
for electroweakly-interacting sparticles via trilep¬ 
tons and dileptons, and for light stops, separately 
using dedicated algorithms validated using the 
Atom and Scorpion 20 codes. In all cases we 
consider the latest SUSY searches from ATLAS 
and CMS that are based on the full Run 1 data 
set, as detailed later in the paper. We perform 
extensive validations of the applications of these 
searches to the pMSSMlO, so as to ensure that 
we make an accurate and comprehensive set of 
implementations of the experimental constraints 
on the model. 

More information about the scan of the 
pMSSMlO parameter space using the MultiNest 
technique, as well as details about our implemen¬ 
tations of the LHC searches, are provided in Sec¬ 
tion [21 Section [3] discusses the results of the 
pMSSMlO analysis, including the best-fit point 
and other benchmark points with low sparticle 
masses that could serve to focus analyses at Run 2 
of the LHC. Section [3] discusses the extent to 
which the preferred ranges of pMSSMlO param¬ 
eters permit renormalization-group extrapolation 
to GUT scales. Section analyses the prospects 
for discovering SUSY in future runs of the LHC, 
Section analyses the prospects for discovering 
SUSY at possible future e+e“ colliders, and our 
conclusions are summarised in Section Ul 


2. Method 

We describe in this Section how we perform a 
global fit of the pMSSMlO taking into account 
constraints from direct searches for SUSY par¬ 
ticles, the Higgs boson mass and rate measure¬ 
ments, SUSY Higgs exclusion bounds, precision 
electroweak observables, H-physics observables, 
and astrophysical and cosmological constraints on 
cold dark matter. We describe the scanned pa¬ 
rameters and their ranges, the framework that 
we use to calculate the observables, and the treat¬ 
ment of the various constraints. 


2.1. Parameter Ranges 

As described above we consider a ten- 
dimension subset (pMSSMlO) of the full pMSSM 
parameter space. The selected SUSY parameters 
were listed in Eq. Q , and the ranges of these pa¬ 
rameters that we sample are shown in Table 
We also indicate in the right column of this Ta¬ 
ble how we divide the ranges of most of these 
parameters into segments, as we did previously 
for our analyses of the CMSSM, NUHMl and 


NUHM2 16 21 


The combinations of these segments constitute 
boxes, in which we sample the parameter space 
using the MultiNest package 17 . For each box, 
we choose a prior for which 80% of the sample has 
a flat distribution within the nominal range, and 
20% of the sample is outside the box in normally- 
distributed tails in each variable. In this way, our 
total sample exhibits a smooth overlap between 
boxes, eliminating features associated with box 
boundaries. An initial scan over all mass param¬ 
eters with absolute values < 4000 GeV showed 
that non-trivial behaviour of the global likelihood 
function was restricted to |Mi| < 500 GeV and 
^ 1000 GeV. In order to achieve high resolu¬ 
tion efficiently, we restricted the ranges of these 
parameters to |Mi| < 1000 GeV and Q < mi < 
2000 GeV in the full scan. 


2.2. MasterCode Framework 

We calculate the observables that go into 
the likelihood using the MasterCode frame¬ 
work 16 21-25 , which interfaces various public 

for 
for 


and private codes: SoftSusy 3.3.9 
the spectrum, FeynWZ 
precision observables, FeynHiggs 2.10.0 
sector and {g — 
for 


29 


for the 


SuFla 


30 


Higgs 
Superlso 


31 


the electroweak 

m 

the B- 
[321 for 


physics observables, Micromegas 3.2 
the dark matter relic density, SSARD 
for the spin-independent cross-section 
for 


SDECAY 1.3b 34 


calculating sparticle 


branching ratios, and HiggsSignals 1.3.0 35 


and HiggsBounds 4.2.0 36 for calculating con¬ 
straints on the Higgs sector. The codes are linked 
using the SUSY Les Houches Accord (SLHA) [3^. 
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Parameter 

Range 

Number of 
segments 

Ml 

(-1,1) TeV 

2 

M 2 

(0,4) TeV 

2 

M 3 

{A, A) TeV 

4 

mq 

(0,4) TeV 

2 


(0,4) TeV 

2 

mj 

(0,2) TeV 

1 

Ma 

(0,4) TeV 

2 

A 

(-5,5) TeV 

1 

9- 

(-5,5) TeV 

1 

tan/3 

( 1 , 60) 

1 

Total number of boxes 


128 


Table 1 

Ranges of the pMSSMlO parameters sampled, together with the numbers of segments into which each 
range was divided, and the corresponding number of sample boxes. 


2.3. Electroweak, Flavour, Cosmological 
and Dark Matter Constraints 

For many of these constraints, we follow very 
closely our previous implementations, which were 
summarized recently in Table 1 in [16| . Specif¬ 
ically, we treat all electroweak precision ob¬ 
servables, all i?-physics observables (except for 
BR{Bs ci —t {9 ~ 

sity as Gaussian constraints. The contribution 
from BR(i3s,(i —t pf'pc), combined here in the 
quantity |^, is calculated using the com¬ 
bination of CMS and LHCb results de¬ 
scribed in [^. We incorporate the current world 
average of the branching ratio for BR (6 —>■ sy) 
from [38] combined with the theoretical estimate 
in the SM from 39 , and the recent measure¬ 
ment of the branching ratio for BR(Ru —>■ tv^) 
by the Belle Collaboration combined with 
the SM estimate from 41 . We use the upper 


limit on the spin-independent cross section as a 
function of the lightest neutralino mass rn^o from 
which is slightly stronger than that 
taking into account the 


LUX 42 


from XENONIOO 43 


theoretical uncertainty on as described in 


21 


2.4. Higgs Constraints 

We use the recent combination of ATLAS and 
CMS measurements of the mass of the Higgs bo¬ 


son: Mh = 125.09±0.24 GeV 44 , which we com¬ 


bine with a one-cr uncertainty of 1.5 GeV in the 
FeynHiggs calculation of Mh in the MSSM. 

In addition, we refine substantially our treat¬ 
ment of the Higgs boson constraints, as com¬ 
pared with previous analyses in the MasterCode 
framework. In order to include the ob¬ 
served Higgs signal rates we have incorporated 
which evaluates the con- 


35 


HiggsSignals 
tribution of 77 channels from the Higgs bo¬ 
son searches at the LHC and the Tevatron (see 


Ref. 35 for a complete list of references). A dis¬ 


cussion of the effective number of contributing 
channels is given in Sect. 13.21 below. 

We also take into account the relevant searches 
for heavy neutral MSSM Higgs bosons via the 
H/A 


channels 


45 


46 


We evaluate 


the corresponding contribution using the code 
which includes the latest CMS 
based on ~ 25 fb“^ of data^ These 


HiggsBmmds 
results 


45 


results include a combination of the two possi¬ 
ble production modes, gg —>■ H/A and bb —> 
bbH/A, which is consistently evaluated depend¬ 
ing on the MSSM parameters. Their implementa¬ 
tion in HiggsBounds has been tested against the 
published CMS data, and very good qualitative 


and quantitative agreement had been found 47 


Other Higgs boson searches are not taken into 


^The corresponding ATLAS results [46| have similar sen¬ 
sitivity, but are documented less completely. 
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account, as they turn out to be weaker in the 
pMSSMlO that we study. 

2.5. LHC Constraints on Sparticle Masses 

A comprehensive and accurate application of 
the SUSY searches with the full Run 1 data 
of the LHC to the pMSSMlO parameter space 
is a central part of this paper. As most of 
these searches have been interpreted by AT¬ 
LAS and CMS only in simplified model frame¬ 
works, we have introduced supplementary pro¬ 
cedures in order to apply these searches to the 
complicated sparticle spectrum content of a full 
SUSY model such as the pMSSMlO. For this 
we consider three separate categories of parti¬ 
cle mass constraints that arise from the LHC 
searches: a) generic constraints on coloured spar- 
ticles (gluinos and squarks), b) dedicated con¬ 
straints on electroweakly-interacting gauginos, 
Higgsinos and sleptons, c) dedicated constraints 
on stop production in scenarios with compressed 
spectra. We refer to the combination of all 
these constraints from direct SUSY searches as 
the LHC8 constraint, with sectors labelled as 
LHCScoi, LHCSewk, and LHCSgtop, respectively. 
In the following subsections we provide further 
details about our implementations of these indi¬ 
vidual constraints, discussing in detail the vali¬ 
dations of our procedures and the corresponding 
uncertainties. 

We use two dedicated software frameworks for 
recasting the LHC analyses used in this paper. 
Both frameworks implement the full list of cuts 
of a given experimental search to obtain yields in 
the respective signal regions of the search. These 
signal yields are then confronted with the SM 
background yields and observations in data, as 
reported by the experimental searches. Based on 
these comparisons we construct the standard sta¬ 
tistical estimator CL^ 48 , which is also used by 


of LHC detectors by mapping from the truth- 


the experiments to determine the compatibility of 
their data with a given signal hypothesis. In this 
way it is possible to interpret the various LHC 
searches in any given SUSY model, such as those 
explored in our pMSSMlO scans. 

To recast the ATLAS searches considered in 


this paper we use Atom 19 , which is a Rivet 49 
based framework. Atom models the resolutions 


level particles found for example in PYTHIA 6 50 


event samples to the reconstructed objects, such 
as 5-jets and isolated leptons, according to the re¬ 
ported detector performances. In particular, the 
efficiencies of object reconstruction and the pa¬ 
rameters associated with the momentum smear¬ 
ing are implemented in the form of analytical 
functions or numerical grids. The program has 
already been used in several studies [^, and the 


validation of the code can be found in . 

For the CMS searches we use a private code 


called Scorpion 20 that was already used in 18 


Scorpion obtains signal yields for a number 
of CMS searches based on events generated 
with PYTHIA 6 50 that are passed through the 


DELPHES 3 53 detector simulation package us¬ 
ing an appropriate data card to emulate the re¬ 
sponse of the CMS detector. A significant effort 
was made to validate the modelling of these anal¬ 
yses by comparing the results obtained with the 
published results of the experimental collabora¬ 
tion. For further information on the validation of 


the CMS searches see 18 


The signal yields from Atom and Scorpion are 
confronted with the background yields and obser¬ 
vations obtained from the individual ATLAS and 
CMS searches, and the corresponding CL^ is cal¬ 
culated using the LandS package . We convert 
the calculated CL^ value for a generic spectrum in 
the MSSM into a contribution by interpreting 
it as a p-value for the signal hypothesis assuming 
one degree of freedom. 

2.5.1. LHC constraints on coloured sparti- 
cles 

In the cases of the CMSSM, NUHMl and 


NUHM2, we showed in 16 21 that it was suf¬ 


ficient to extrapolate to other parameter values 
the exclusion contour in the CMSSM (mo,mi/ 2 ) 
plane from the ATLAS search for jets-h^r 
that was given for specified values of tan/? and 
Aq. We showed that the ATLAS exclusion is, 
to good approximation, independent of tan /? 


and Aq 16 23 ^ and, for the applications to 


the NUHMl and the NUHM2, we checked that 
these limits in the (mo,TOi/ 2 ) plane were inde¬ 
pendent of the degrees of non-universality of the 
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soft SUSY-breaking contributions to the Higgs 
masses, within the intrinsic sampling uncertain¬ 
ties. 

In the case of the pMSSMlO, however, the im¬ 
plementation of the direct searches for coloured 
sparticles is less straightforward. It is computa¬ 
tionally impossible to apply all the LHC search 
constraints individually to each of the ~ 1.2 x 10 ® 
parameter choices in our sample. For example, 
PYTHIA 6 and DELPHES 3 take several minutes for 
the generation of 10,000 events followed by detec¬ 
tor simulation, which is required to determine the 
signal acceptance and CLg of each point sampled 
in the parameter space. Instead, we follow an 
approach outlined in 18 , which constructs uni¬ 


lish lower limits on the gluino mass, rrig, 
third-generation squark mass, niq^, 


versal mass limits on coloured sparticles by com¬ 
bining an inclusive set of jets -I- X +1^t searches, 
as we now describe. 

As was shown in [Ts], it is possible to estab- 

and the 
that are in¬ 
dependent of the details of the underlying spec¬ 
trum, within the intrinsic sampling uncertainties, 
by combining a suitable set of inclusive SUSY 
searches. In this approach the limits only depend 
on rrig, rriq^ and the mass of the lightest sparticle 
m^o. The essence of the idea is that strongly- 
interacting sparticles decay through a variety of 
different cascade channels, whose relative proba¬ 
bilities depend on other model parameters. How¬ 
ever, if one combines a sufficiently complete set of 
channels of the form jets -p X +1^t, one will cap¬ 
ture essentially all the relevant decay channels. 

In order to apply this idea to the pMSSMlO pa¬ 
rameter space, we have to extend this approach 
to include also the generic first- and second- gen¬ 
eration squark mass, m^, as a free parameter. We 
then construct a ‘universal’ function that de¬ 
pends only on m^o, rrig, rriq, and as de¬ 

tailed below. This function defines our imple¬ 
mentation of this LHCScoi constraint. There are 
two caveats to this approach. One is that the re¬ 
gion of parameter space where — m^o is small, 
which is the object of dedicated searches, requires 
special attention. The other is that searches for 
electroweakly-produced sparticles (sleptons, neu- 
tralinos and charginos) fall outside the scope of 
the LHCScoi constraint. We have developed ded¬ 


icated approaches to establish accurate LHC lim¬ 
its for the special cases of electroweakly-produced 
sparticles and the compressed-stop scenario with 


2.5.2 


— m^ o < mt, as described in Sections 
and 2.5.3| respectively. 

In order to construct y® as a function 

we first 


of m^o , m 


X 

and TOq 3 , 


g, rriq, and 771 ^ 3 , we nrst gener¬ 
ate a sample of points on a 1 -I- 3 dimen¬ 
sional grid, which we use for linear interpo¬ 
lation. We construct this grid starting from 
= {10, no, ..., 610} GeV. 

we select 


m^o, 


{m^o 


values of m^o 

/V 2. 

For each of these values of 
the following values of rrig and niq. 

40, m^o + 140, ..., 1750, 2500, 5000} GeV, 
whereas 777^3 takes values { 777^^0 + 80, m^o + 
180, ..., 1290, 2500, 5000} GeV, where \he 
dots indicate steps of 100 GeV, so that the total 
number of points in the grid is 25,564. The choice 
for this grid is motivated by the need for a fine 
granularity at low masses, while also capturing 
the parameter behaviours at higher masses. 

We associate a SUSY spectrum to each point 
on the grid, by setting the first- and second- 
generation squark masses equal to rriq, and the 
third-generation squark masses equal to 777 ^ 3 . For 
each SUSY spectrum we generate coloured spar¬ 
ticle production events using PYTHIA 6 50 and 


pass them through the DELPHES 3 53 detector 


simulation code using a detector card that emu¬ 
lates the GMS detector response. We then pass 
the resulting events through Scorpion |^, which 
emulates the monojet, MT2, single-lepton, same- 
and opposite-sign dilepton (SS and OS) and 3- 
lepton GMS searches 56, ^ , to estimate the 


numbers of signal events in each of the signal re¬ 
gions. After this we calculate the GLg using the 
Lands package 54 , by combining all signal re¬ 


gions from these searches. If searches have over¬ 
lapping signal regions, we take the strongest ex¬ 
pected limit, as is the case for the GMS monojet 
and single-lepton searches. 

In Fig. [^we show a three-dimensional overview 
and a pair of two-dimensional slices through 
this grid. The top panel shows the full three- 
dimensional grid for rn^o = 310 GeV and illus¬ 
trates the fine and coarse granularity of the grid 
at low and high values of rrig rriq, and 777 ^ 3 , re¬ 
spectively. The lower left panel shows the two- 
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Figure 1. Illustration of the grid in ni^o, nig, rriq, and 771^3 on which x^(Scorpion) is evaluated in order to 
construct LHCScoi- The upper panel shows the three-dimensional grid for m^o = 310 GeV, the lower left 
panel shows a two-dimensional slice through the grid, and the lower right panel is another two-dimensional 
slice that illustrates the dependence on ruj^o, see the text. 


dimensional slice for the same neutralino mass 
and mq^ = 2500 GeV, highlighting that there 
is only a small, though non-negligible, depen¬ 
dence of the function on mq for values of 
TOg > 2500 GeV. The lower right panel shows 
the function as a function of mq„ and m^o , for 
fixed mq = 2500 GeV and nig = 2500 GeV, illus¬ 
trating that for different values of different 
grids are defined in mg, m^, and 

In order to apply the LHGScoi constraint to 


a generic pMSSMlO spectrum, we calculate to, 
( 771 ^ 3 ) as the cross-section-weighted average of 
the first- and second- (third-)generation squark 
masses, to ensure that the LHCScoi constraint re¬ 
flects the actual production cross-sections. This is 
especially relevant for the third-generation squark 
masses, as they generally have large splittings. 
The contribution for LHCScoi is obtained by 
linear interpolation of the x^ values on the 1-1-3- 
dimensional grid. There is one special case when 


































9 


mf_^ —m^o < mt'. here the standard searches listed 
above are less sensitive, and the universality of 
the limits is expected to break down. In this case, 
we calculate rriq^ assuming zero cross-section for 
the lighter stop, and consider separately the im¬ 
pacts of dedicated stop searches in this region, as 
described in Section 12.5.31 

In order to validate the LHCScoi constraint and 
to gauge quantitatively its uncertainty, we have 
performed a number of studies and tests. First, 
we randomly selected 1000 model points from our 
sample where at least one of the sparticle masses 
is low enough to have been within the reach of 
LHC Run 1 (m^o < 600 GeV and either rrig < 
1500 GeV, rriq <'l600 GeV or rriq^ < 900 GeV) 
and < 10 relative to the global minimum. 
For these points we compare the values in¬ 
terpolated from the look-up table (y^ (LHCScoi)) 
with the y^ obtained by running the full chain of 
event generation, detector simulation and anal¬ 
yses (y^(Scorpion)). The left panel of Fig. 
shows a histogram of the differences for the 1000 
randomly-selected points. As indicated in the leg¬ 
end of this figure, the standard deviation on this 
distribution is a ^2 = 1.8. 

The right panel of Fig. [^shows a scatter plot in 
the (y^(Scorpion), y^(LHC8coi)) plane of the y^ 
values obtained from the two approaches. They 
would agree perfectly along the diagonal where 
y^(Scorpion) = y^(LHC8coi), and the lighter- 
and darker-shaded blue strips are the ±la ^2 and 
±2 ct^2 bands around this diagonal. The verti¬ 
cal and horizontal dashed lines in this plot corre¬ 
spond to the 95% CLg in each approach. For the 
majority of points, the interpolation and the full 
analysis agree whether the point is excluded at 
the 95% CLg, or not, and most of the remaining 
points lie within ± 2 a^ 2 . 

We then assess how the uncertainty cr ^2 in our 
implementation of the LHC8coi constraint trans¬ 
lates into uncertainties in sparticle mass limits: 
see the upper left panel of Fig. For this esti¬ 
mate, we bin the 1000 points of the first test, and 
calculate the standard deviation, a^ 2 , for points 
with y^ (LHCScoi) < 1, 1 < (LHCScoi) A 4 and 

^The other panels of Fig.|^show the corresponding uncer¬ 
tainties in our treatments of the LHCSbwk and LHCSstop 
constraints, which are discussed later. 


y^ (LHCScoi) > 4. We then apply the LHCScoi 
constraint in three ways: with the nominal imple¬ 
mentation, and shifting the y^ (LHCScoi) penalty 
up and down according to these binned standard 
deviations. The results are shown in the upper 
left panel of Fig. as solid and dotted red (blue) 
contours in the {mg,mq) plane corresponding to 
the nominal and up- and down-shifted cases for 
the 68 (95)% CL, respectively A dedicated 
study of points within the 68% and 95% CL re¬ 
gions confirms that our implementation of the 
LHCScoi constraint is valid within these uncer¬ 
tainties, and our estimate of y^ at the best-fit 
point differs from the Scorpion evaluation by less 

than one El 

We conclude that the uncertainty cr ^2 in our 
estimate y^ (Scorpion) is generally reliable, and 
translates into an uncertainty of 0(50 GeV) in 
the limits on the gluino and squark masses, which 
is fully sufficient for the purpose of our studies. 


2.5.2. LHC constraints on electroweak 
gauginos, Higgsinos and sleptons 

Unlike the searches for coloured sparti- 
cles, where we were able to construct a 
computationally-efficient, approximately univer¬ 
sal limit, the LHG constraints on electroweakly- 
produced sparticles vary strongly in sensitivity, 
depending on the mass hierarchy of sparticles 
and their corresponding decay modes and fi¬ 
nal states. For example, searches in the three- 
lepton plus missing energy channel constrain the 
chargino and neutralino masses up to m~± = 

Xl 


m^o < 700 GeV for m^o < 300 GeV, if yf and ya 




decay exclusively into on-shell sleptons 58 


59 


whereas a much weaker limit, m-± = m^o ^ 450 

’ Xl X2 ~ 


GeV for m^o < 100 GeV, was found in an anal- 
ysis of the two-lepton plus missing energy chan- 
nel 59 60 , assuming that the xt X2 decay 
exclusively into the Xi association with W and 
respectively, and not taking into account the 


^Here and in subsequent analogous parameter planes, we 
treat the = 2.30 and 5.99 contours as proxies for the 
68% and 95% CL contours. 

“^Our LHCSewk and LHCSstop analyses described later 
also differ by less than one from the corresponding 
Scorpion/Atom evaluations. This is also true for the 
benchmark points introduced later. 













10 



id {LHC8^„i)~x^ (Scorpion) 



Figure 2. Left panel: Histogram of the differences between the values of the likelihood function 
(Scorpion) evaluated using individual LHCScoi searches for 1000 randomly-selected points and the 
estimate (LHCScoi) obtained by interpolation from a look-up table as described in the text. Right panel: 
Scatter plot in the (x^(Scorpion), x^(LHC8coi)) plane of the x^ values obtained from the two approaches; 
the vertical and horizontal dashed lines in this plot correspond to the 95% CLg in each approach. 


decay X 2 Xi^ 61,62 . The same two-lepton 


analyses constrain slepton pair production, lead¬ 
ing to the limits ^ ^ 270 (200) GeV for 

m^o < 100 (50) GeV f^j^. Therefore, the uni¬ 
versal limit approach tliat we use to combine and 
characterise searches for coloured sparticles is in¬ 
applicable to searches for electroweakly-produced 
sparticles, and we use an alternative method. 

For model points where the production of 
electroweakly-produced sparticles provides a non¬ 
trivial constraint, they must be much lighter than 
the coloured sparticles, since otherwise the much 
higher rates of production of coloured sparticles 
would already exclude the model points. There¬ 
fore, in the region of interest, there can be only 
a few particles lighter than the electroweakly- 
produced sparticles, implying that one can use a 
combination of a few simplified models (SMS) to 
approximate the sensitivities of the LHG searches 
for the production of these sparticles. Depend¬ 
ing on the decay mode and final state, we select 
ATLAS and/or GMS limits derived from relevant 
simplified models to calculate the contributions 
of these searches to our global x^ function. For 


the LHC searches that constrain electroweakly- 
produced gauginos, Higgsinos and sleptons, to 
a good approximation all relevant x^ contribu¬ 
tions can be extracted from simplified chargino- 
neutralino and simplified smuon and selectron 
models. 

For each simplified model limit we construct 
a function Xsms depends on the two rele¬ 
vant masses: (m~± ~ m^o,m^o) for the simpli- 

^ Xi X2 ’ X.\' ^ 

fied chargino-neutralino model and {m^, m^o) for 
the simplified slepton (£ = e,fi) model. We as¬ 
sume that Xsms = 15 in the bulk of the region 
excluded in the simplified model, and that this 
X^ penalty vanishes exponentially when crossing 
the boundary to the allowed region, with the gen¬ 
eral form 


Xsms 


min 

l,r 


15- B- 


1 


( 2 ) 


where the subscripts I, r indicate the simplified 
model exclusion contour to the left and right (in 
the horizontal direction, i.e., ~ m^o or mf) 

of the point on the contour with the largest value 
of ruyo, B is the branching ratio of the decay in 
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Figure 3. Impacts of the ± 1 a uncertainties in our implementations of the LHCScob LHCSewk 
LHCSstop constraints on the 68 and 95% CL regions (indicated by the red and blue contours) in the cor¬ 
responding relevant mass planes: {mq,mg) (upper left panel), {m^±,m^o) (upper right panel), 

(lower left panel), and (rnf_^,m^o) (lower right panel). In each case, the dot-dashed and dashed contours 
are obtained by shifting the respective penalty up and down by one standard deviation , as discussed 
in the text. The filled green stars correspond to the nominal best-fit point and the open stars (shown if 
not overlapping) to those which were obtained from shifting the up or down with cr^ 2 . We note that 
in the lower right panel the best-fit points lie outside the displayed parameter range. 


question (as calculated with SDECAY [^), d is 
the closest distance in GeV to the contour, and 
H and a control the precise fall-off of the func¬ 
tion, so as to mimic the experimental uncertainty 
bands, and are functions of m.^o. We note that if 
one sets fi = —a then xiMs('^ = 0) « 4, so that 
the exclusion on the contour corresponds approx¬ 
imately to the 95% CLg. Finally, to avoid an un- 
physically slow fall-off outside the 95% CL^ limit 
we set cr = 50 GeV and adjust d accordingly if 


cr > 50 GeV and d—fj, > a (and hence Xsms ^ 

In order to illustrate Eq. ([^, we display in 
Fig-ElxiMs/^ the xfx 2 decay via sleptons. 
In the left panel XsMs/dd is shown for a fixed value 
of m^o = 300 GeV where the green (blue) line 
corresponds to di, p,i,ai, {d^-, whereas ver¬ 

tical dashed lines indicate the position of the con¬ 
tour. The right panel shows the same Xsms/-® 
colour) as a function of w-± ~ m^o and m^o, and 
the 95% GLg exclusion contour found in Fig. 7(a) 
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of (blue line). Note that we apply no con¬ 
straint for m^o > 380 GeV, the highest value on 
the blue experimental contour. 

In order to establish LHCSewk we tuned the p, 
and cr parameters for each simplified model to re¬ 
produce best the values that we obtained using 
Atom for a representative set of model points from 
our sample. Table [^summarises the implementa¬ 
tions of the simplified model exclusion limits that 
contribute to LHCSewk- Note that, as described 
above, the large value of ar = 300 GeV for the 
limit from xfx 2 production and decay via WZ 
is replaced by setting Uj. = 50 GeV and adjust¬ 
ing dr accordingly when dr — > o'r (and hence 

XsMS ^ Also, we had to produce our own con¬ 
tour for the direct production of right- and left- 
handed sleptons (selectrons and smuons), corre¬ 
sponding to their production cross-sections. Note 
that this simplified model contour is also applied 
when left-handed sleptons decay via X 2 xf ■ 

In order to validate our method and to deter¬ 
mine quantitatively its uncertainty, we compare 
the contributions to the global x^ function calcu¬ 
lated with this LHCSewk limit approach, 

(LHCSewk) = ^ Xsms > (3) 

SMS 

to results from a full recast of all the above- 
listed searches as implemented in Atom. In this 
recast the full analysis is simulated, so that it 
is possible to determine for any arbitrary SUSY 
spectrum the CL^ value (and hence the corre¬ 
sponding x^) with which a given search penal¬ 
izes the SUSY spectrum. We obtain a set of 
1000 model points from our sample by binning the 
(m^o « I pl^'ae in 100 x 100 bins, select¬ 
ing one point randomly per bin, and then take a 
random subset of 1000 of these points. This pro¬ 
cedure was employed to ensure a representative 
set of the decay modes in our sample. 

Fig. [^displays scatter plots in the (to-±,to^o) 
plane of the contributions to the global func¬ 
tion for these 1000 model points as calculated us¬ 
ing the LHCSewk method (y^(LHCS ewk)) (left 
panel) and the Atom code y^(Atom) (right panel), 
with the indicated colour code in each plot. The 
immediate visual impression is that the colours in 


the two scatter plots are generally quite similar, 
indicating that the two procedures deliver similar 
y^ contributions overall. A closer inspection of 
the plots reveals similar bands of low-y^ points 
with small ni-± — m^o in a chargino coannihila¬ 
tion strip region, while elsewhere we see similar 
disfavouring of points with low m^o ^ 150 GeV 
and larger m~±. However, even within this band 
we see a sparse set of points with relatively low 
y^ that appear similarly in both the LHCSewk 
analysis based on simplified models and the Atom 
implementation of the full searches. These are 
mainly due to the decay X 2 Xi^, thus weaken¬ 
ing the stronger y^ —y^Z-based limit. 

For a more quantitative comparison of our 
LHCSewk method and Atom we turn to Fig. [^ 
We see in the left panel that the difference 
between y^ (LHCSewk) and y^(Atom) is rela¬ 
tively small, with an r.m.s. difference a ^2 = 

2.31. The correlation between y^(LHCS ewk) and 

y^ (Atom) is visible in the scatter plot in the right 
panel of Fig. We see that most points are ei¬ 
ther excluded with Ay^ > 4 in both analyses, or 
allowed with Ay^ < 4 in both cases. Last but 
not least, there are relatively few ‘off-diagonal’ 
points with large Ay^, which form the small non- 
Gaussian tail of the y^ (LHCSewk) — X^(Atom) 
distribution seen in the left panel of Fig. [^ 

To quantify the impact of this uncertainty on 
our analysis, we follow the same procedure as for 
our limits on coloured sparticles, and translate 
the cr ^2 (binned analogously) into a ± 1 cr band 
for our 68% and 95% CL contours in the impor¬ 
tant {m^± ,m-^o) and (m^^,mjjo) planes. As can 
be seen in the upper right and lower left panels of 
Fig.i the uncertainty associated with LHCSewk 
is in general small in the 68% CL region of our 
fit, although it is larger at the 95% CL level in 
the (TO-±,mj^o) plane. The effects on the best-fit 
point of these upward and downward shifts in the 
y^ treatment are shown in these panels as open 
green stars. The downward shift has very little ef¬ 
fect, and is essentially invisible in the (to-± , m^o) 
plane. The upward shift increases the best-fit 
values of m^o and m~± while reducing that of 

Xl Xl ^ 

though the variations are contained well 
within the 68% CL region, clearly indicating that 
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Figure 4. Illustration of Xsms/^^ defined in Eq. for xfx 2 production and decay via sleptons. 
In the left panel xIms/-^ shown for a fixed value of m^o = 300 GeV, where the green (blue) line 
corresponds to di,fj,i,ai, (dr, fJ.rjO'r) and vertical dashed lines indicate the position of the contour. The 
right panel shows the same Xsms/^ colour) as a function of m^± ~ m^o and m^o, and the 95% CL^ 
exclusion contour found in Fig. 7(a) of \5^ (blue line). 


Simplified Model 

Limit 

[p-uai] [GeV] 

[GeV] 

X 1 X 2 via i 

Fig. 7(a) in 

58 


(-5, 5) 

(-40, 40) 

xfx2 via WZ 

Fig. 7(b) in 

58 


(-20, 20) 

(-300, 300) 


Generated using Atom 

(-20, 10) 

(-40, 30) 


Table 2 

The simplified model limits used to constrain electroweak gauginos, Higgsinos and sleptons. 


the corresponding uncertainties do not impact the 
overall conclusions. 


2.5.3. LHC constraints on compressed stop 
spectra 

In their searches for stop production, ATLAS 
and CMS have placed special emphasis on com¬ 
pressed spectra, which pose particular challenges 
for LHC searches. Whilst limits on stop produc¬ 
tion in the region where — m^o > mt are 
fully included in the LHCScoi limits described 
in Section |2.5.1[ a dedicated treatment of the 
compressed-spectrum region — m^o < mt is 
required in order to include properly all the rel¬ 
evant collider limits. In this region we calculate 
the contribution of stop searches to the global x^ 
in a similar way as for the for electroweakly pro¬ 
duced sparticles described in Section 2.5.2. We 


refer to this dedicated limit-setting procedure as 

LHC8,top. 

We show in Fig. [^a colour-coded scatter plot 
in the (m(^,TO^o) plane of the ti decay modes 
with branching ratios > 50% for 1000 randomly- 
selected pMSSMlO points in the region of inter¬ 
est. We see that the ii —)■ bx^ mode (shown in 
light green) dominates for the majority of points, 
and that this decay can be important through¬ 
out the parameter region displayed. We also find 
that, when this is the dominant stop decay mode, 
in most cases the xf and xj are almost mass de¬ 
generate. To constrain the final states with this 
decay mode we implement the simplified model 
limit presented in Fig. 6 of the ATLAS di-bottom 
analysis 63 , where w-± — m^o = 5 GeV is as¬ 
sumed, applying this for the model points with 
m-± — TOvO < 30 GeV. 

Xi Xi 
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Figure 5. Scatter plots in the {m^±,m^o) plane of the contributions to the global functions from 
the electroweakly-interacting sparticle constraints for 1000 randomly-selected points accessible to LHC 
searches, as calculated using the LHCSewk method based on simplified model searches (x^(LHC8ewk), 
left panel) and the Atom code right panel). 



0 (LHCSewk)(A tom) 



Figure 6. Left panel: Histogram of the differences between the values of the contributions of the 
electroweakly-interacting sparticle constraints to the global likelihood function (LHCSewk) evaluated 
using simplified model searches for the 1000 randomly-selected points and the estimate x^(Atom) obtained 
using the Atom code. Right panel: Scatter plot in the (x^(Atom), x^(LHC8ewk)) plane of the values 
obtained from the two approaches; the vertical and horizontal dashed lines in this plot correspond to the 
95% CLs in each approach. 


(Atom) 
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Figure 7. Scatter plot in the plane of the ti decay modes with branching ratios > 50% for 

1000 randomly-selected points with — rn^o < mt- 


If — m^o > M^r + TOb, the 3-body ti —>■ 
hWffl mode can dominate stop decay. The points 
for which this mode is dominant are shown by 
purple dots in Fig. For this decay mode we 
implement the simplified model limit presented 
for Mw +rnb < — m^o < mt in Fig. 15 of the 


ATLAS single-lepton analysis 64 

In the — m^o < M\y region, the 


decays ii —)• cxi (red dots in Fig. and ii —)• 
bff'Xi (grey dots) can be the dominant stop de¬ 
cay modes. The ii —)• bi'rTi mode (green dots) 
may also dominate stop decay in this region, as 
well as in the m^_^ — ^ Mw mf, region, as 

can also be seen in Fig. 

Due to the variety of different stop decay modes 
that are relevant in this compressed region, we 
cannot use only the limits from simplified models 
provided by the experiments, as they do not cover 
all relevant decay chains and assume branching 
ratios of 100%. However, these missing, in part 
rather complex, decay chains can effectively be 
constrained by hadronic inclusive searches such 
as those we have already used for our LHCScoi 


limits. In particular, the CMS hadronic mT 2 
search 56 has rather high sensitivity for these 
decay chains, as the kinematic phase space cov¬ 
ered by the search makes no special assumptions 
on the final state, other than it having a purely 
hadronic signature. 

Based on these inclusive searches, we derive 
limits for simplified models for ti —)■ ex)* and 
ti —>■ biyrfi decays. For the ti —> sim¬ 

plified model we assume — m^o < 40 GeV 
when creating the limit in the (m^^, rn^o ) plane. 
We do not implement a simplified model limit for 
ii —t &//'xi because this decay mode has negligi¬ 
ble impact on our study, as can be seen in Fig. 
Using these simplified model limits, we constrain 
the stop decay modes following a procedure very 
similar to what we used for LHCSewK; using an 
interpolating function of the form to mimic 
the uncertainty (yellow) band in, e.g., Fig. 6c 


63 . We summarise our implementation of the 


simplified model limits in Table 13 When estab¬ 
lishing these limits we use values of the param¬ 
eters and apr that depend on m^o. When- 
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ever multiple values of these parameters are given 
for different values of m^o, the parameters for in¬ 
termediate values of rn^o are obtained by linear 
interpolation, and taken as constants elsewhere. 

As for our LHCScoi and LHCSewk limit im¬ 
plementations, it is also important to deter¬ 
mine accurately the uncertainty in the dedicated 
limit procedure for the compressed stop region. 
Note that in the compressed region not only the 
constraints from LHCSstop but also those from 
LHCSewk play a role. Therefore we first assess 
the qualitative agreement between (LHCSstop) 
and the “true” x^(Atom and Scorpion) as cal¬ 
culated using the Scorpion and Atom codes, 
for points with x^(LHCSewk) < 2. Fig. 
compares scatter plots in the {mi _^, m^o) plane 
of (LHCSewk) + x^(LHCSstop) (left panel) 
and x^(Atom and Scorpion) (right panel). The 
colour code used is indicated on the right-hand 
sides of the panels, and we see that the patterns 
of colours in the two scatter plots are qualitatively 
similar. This is remarkable, given the interplay of 
so many different decay chains. 

More quantitative comparisons of the con¬ 
tributions to the global x^ function calculated 
on the basis of the simplified model searches 
for stops and electroweakly produced sparti- 
cles (x^ (LHCSewk) + x^ (LHCSstop)) with re¬ 
sults from Scorpion and Atom for these 1000 
randomly-selected pMSSMlO points (x^(true)) 
are shown in Fig. The left panel shows a his¬ 
togram of the difference between x^ (LHCSewk ) + 
X^(LHCSstop) and x^(ti'ue), showing that it is rel¬ 
atively small, with an r.m.s. difference a ^2 = 3.15. 
The right panel of Fig. displays a scatter plot 
in the (x^(Atom and Scorpion), x^(LHCS ewk) + 
X^(LHCSstop)) plane. We see that points that are 
(dis)favoured at the 95% CLg level in the simpli¬ 
fied approach are, in general, also (dis)favoured 
at the 95% CL^ level in the more sophisticated 
approach based on Scorpion and Atom. 

To determine quantitatively the effect of the 
uncertainty in the LHCSstop procedure, we trans¬ 
late the impact of the above-mentioned a ^2 = 
3.15 uncertainty into the (m(^,m^o) plane in the 
lower right panel of Fig. This shows the im¬ 
pacts of ± 1 a ^2 variations on our 6S% and 95% 


contours in this plane, which is rather small ex¬ 
cept for small values of and m^o. 

Based on this study, we conclude that the 
computationally-manageable simplified approach 
LHCSstop is sufficiently reliable for our physics 
purposes. Specifically, we note that there are 
points with low rrif^ that survive the full LHC 
constraints with relatively low x^- 

3. Results 
3.1. Mass Planes 

Fig.[T^ displays the two-dimensional profile 
likelihood functions in planes of (from top left to 
bottom right) the masses of the gluino, the first- 
and second-generation squarks, the lighter stop 
and sbottom squarks, the lighter chargino and the 
lighter stau, each versus the lightest neutralino 
mass m^o. In each panel the solid (dashed) 
red/blue contours denote the 68%/95% CL con¬ 
tours for the case where we do (not) apply any 
LHC constraints, respectively!^ The green filled 
and empty stars indicate the corresponding best- 
fit points. In the cases of the gluino and squarks, 
the filled stars lie beyond the displayed parts of 
the corresponding planes, and their locations are 
indicated by arrows. In these cases the likelihood 
function varies little as a function of the coloured 
sparticle mass. 

On the other hand, we find that in general 
"ly? ^ GeV at the ^ 68% CL, increasing to 
^ 500 GeV at the ~ 95% CL. This and the pref¬ 
erence for low stau masses (< 700 GeV at the 
- 68% CL, < 1000 GeV at the ~ 95% CL) are 
reflections of the fulfilment of the {g — 2)^ con¬ 
straint in the pMSSMlO, cf.. Fig, [l^ below, and 
(in the latter case) the restriction to a common 
slepton mass for all three generations. 

We can distinguish two ranges of that are 
allowed at the 95% CL: a narrow band where 
"ly? ^ ^ broader region at larger 

m^o that also includes regions favoured at the 
68% CL. In the low-m^o region, before apply¬ 
ing the LHC8 constraints the smuon, selectron 
and stau could have been relatively light, and t- 
channel sfermion exchange could bring the relic 
density into the range allowed by cosmology. 


^However, the LEP SUSY constraints [67| are applied. 
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Decay 

Limit 

m^o [GeV] 

[GeV] 

iHr,crr) [GeV] 

Condition/Remark 

i bxf 

Fig. 6(c) in 

63 


210 

300 

(10, 20) 
(-250, 200) 

(-50, 50) 
(-200, 200) 

m~± — m^o < 30 GeV 

Xi Xi 

i —)• 6ITxi 

Fig. 15 in 

64 


100 

150 

(-20,50) 

(-50, 50) 

(-70, 50) 
(-100,50) 

M\y < — m^o < mt 

t bvTi 

Generated using 
Scorpion 

- 

(-50, 50) 

(-20, 50) 

Based on 56 , assuming 
— m^o < 40 GeV 

i cx\ 

Generated using 
Scorpion 

- 

(-20, 20) 

(-20, 20) 

Based on 56 


Table 3 


The simplified model limits used to constrain scenarios with compressed stop spectra. When establishing 
these limits we use values of fii^r o,nd ai^r in Eq. Q) that in some cases depend on m.^o. Whenever multiple 
values of these parameters are specified for different values ofm^o, the parameters for intermediate values 
of m^o are obtained by linear interpolation, and taken as constants elsewhere. 
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Figure 8. Scatter plots in the plane of the contributions to the global functions from the 

ATLAS mono-jet 16^ and single-lepton [66f searches for 1000 randomly-selected points in the regions of 
interest. The left panel shows calculations using simplified model searches (LHCSewk and the right 
panel shows results from the Scorpion and Atom codes (x^(true)j. 


However, after applying the LHC8 constraints 
only the Z- and h-funnels are allowed in this 
region. In the region where m^o > 80 GeV, 
before implementing the LHC8 constraints stau 
coannihilation and t-channel sfermion exchange 
were both possible. However, after applying the 
LHC8 constraints the dominant processes con¬ 
trolling the dark matter density are X? — X® ~ xf 


coannihilations, with the LSP having mainly a 
Bino composition. 

The two top panels of Fig. [T^ display clearly 
the direct impacts of the LHC8 constraints, which 
are visible in the displacements to larger masses 
of the 68% and 95% CL contours, as can be seen 
from the comparison of the solid and dashed lines. 
On the other hand, the pictures in the two middle 
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(LHCSewk) +X^ (LHCSgtop)—(Atom and Scorpion) x^ (Atom and Scorpion) 

Figure 9. Left panel: Histogram of the difference between the values of the contributions of the stop 
constraints to the global likelihood function x^(LHC8ewk) + X^(LHC8stop) evaluated using simplified 
model searches for 1000 randomly-selected points and the estimates of found using Scorpion and 
Atom. Right panel: Scatter plot in the (x^(true), x^(LHC8ewk) + X^(LHC8stop)) plane of the values 
obtained from the two approaches; the vertical and horizontal dashed lines in these plots correspond to 
the 95% CLs in each approach. 


panels are more complex. There are intermediate 
values of that are disfavoured by the LHC8 
constraints, but there are regions with low values 
of that are allowed by the LHC8 constraints 
at the 95% CL, and even some points with 
and that are favoured at the 68% CL, though 
these are not prominent. In the case of the lighter 
sbottom, the LHC8 constraints disfavour the re¬ 
gion where both and have small values. 
However, a small value of is still allowed at 
the - 95% CL if m^o > 300 GeV to 450 GeV, 
where some points are favoured at the 68% CL. 

Finally, the bottom two panels of Fig. [T0| show 
the impacts of the LHC8 constraints on the 
chargino and stau masses. The main impact on 
the chargino mass is to disfavour most values ex¬ 
cept some where m~± — m^o is small. This is an 

^ Xi Xi 

indirect effect of the LHC8 constraints, with the 
coannihilation of the dark matter particle with 
the lighter chargino playing an important role in 
bringing the dark matter density into the allowed 
range. This compression of the spectrum can be 
attributed to the LHC8ewk limits on direct pro¬ 
duction of light sleptons, and to a lesser extent 


on charginos decaying via sleptons. These con¬ 
straints on light sleptons disfavour the t-channel 
sfermion exchange and stau coannihilation re¬ 
gions. The latter is a consequence of our choice 
of a single mass parameter for the masses of all 
the scalar leptons (see also Sect. [7|. In the case of 
the lighter stau, we see in the bottom right panel 
of Fig. [^a triangular region that is favoured at 
the ^ 68% CL, which is somewhat reduced and 
shifted towards higher mass values by the LHC8 
constraints. 


3.2. The Best-Fit Point 

We now discuss the characteristics of the best- 
fit point, whose parameters are listed in Table HI 
together with the parameters of several bench¬ 
mark points that are discussed below. The best- 
fit spectrum is shown in Fig. El and its SLHA 
file 37 can be downloaded from the MasterCode 
website [2^. We note first the near-degeneracy 


between the Xi , X 2 xf ; which is a general fea¬ 
ture of our 68% CL region that occurs in order to 
bring the cold dark matter density into the range 
allowed by cosmology: see the bottom left panel 
of Fig.fT^ Correspondingly, we see in Tabled 

















* - pMSSMlO w LHC8: best fit, Icr, 2ct ★ - pMSSMlO w LHC8: best fit, Icr, 2a 

* ___ - pMSSMlO w/o LHC8: best fit, Icr, 2cr * - - pMSSMlO w/o LHC8: best fit, Icr, 2cr 




★ - - pMSSMlO w LHC8: best fit, Icr, 2cr 

* - - pMSSMlO w/o LHC8: best fit, Icr, 2cr 


★ - - pMSSMlO w LHC8: best fit, Icr, 2cr 

* - - pMSSMlO w/o LHC8: best fit, Icr, 2cr 




- pMSSMlO w LHC8: best fit, Icr, 2cr 

- pMSSMlO w/o LHC8: best fit, Icr, 2cr 


— pMSSMlO w LHC8: best fit, Icr, 2cr 

- - pMSSMlO w/o LHC8: best fit, Icr, 2cr 




Figure 10. The two-dimensional profile likelihood functions for (top left to bottom right) the masses of the 
gluino, the first- and second-generation squarks, the lighter stop and sbottom squarks, the lighter chargino 
and the lighter stau, each versus the lightest neutralino mass In each panel the solid (dashed) 

red/blue contours denote the Ax^ = 2.30/5.99 level contours for the case where we do (not) apply the 
LHC8 constraints, respectively. The green filled and empty stars indicate the corresponding best-fit points. 
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Parameter 

Best-Fit 

Low m; 

^1 

Low mq 

Low mg 

Low all 

Ml 

170 GeV 

300 GeV 

210 GeV 

190 GeV 

-120 GeV 

M2 

170 GeV 

310 GeV 

220 GeV 

200 GeV 

160 GeV 

M3 

2600 GeV 

1660 GeV 

3730 GeV 

- 1070 GeV 

1700 GeV 

mq 

2880 GeV 

3700 GeV 

1530 GeV 

2430 GeV 

1790 GeV 


4360 GeV 

720 GeV 

1840 GeV 

3780 GeV 

1300 GeV 

mi 

440 GeV 

390 GeV 

430 GeV 

410 GeV 

740 GeV 

Ma 

2070 GeV 

3540 GeV 

2810 GeV 

2990 GeV 

1350 GeV 

A 

790 GeV 

1790 GeV 

2510 GeV 

3000 GeV 

1863 GeV 

A* 

550 GeV 

1350 GeV 

640 GeV 

530 GeV 

190 GeV 

tan/3 

37.6 

37.3 

40.8 

33.9 

35.4 


Table 4 

Parameters of the pMSSMl 0 best-fit point and other comparison benchmark points at low , low mq 
and/or mg. 


that Ml ~ M 2 , though M 3 is very different. The 
overall Xi/xVx? rnass scale is bounded from be¬ 
low by the LEP and LHCSewk constraints, and 
from above by (5 — 2 )^, especially at the 68 % CL. 
We display in Fig. [T| the 95% ( 68 %) CL inter¬ 
vals in our fit for the masses of pMSSMlO parti¬ 
cles as lighter (darker) peach shaded bars, with 
the best-fit values being indicated with blue hor¬ 
izontal lines 12 Turning back to Fig. m we note 
the near-degeneracy between the slepton masses, 
which reflects our assumption of a common in¬ 
put slepton mass at the input scale Msusy that 
would not hold in more general versions of the 
pMSSM. The overall slepton mass scale is below 
1 TeV, as seen in Fig. being bounded from 
above by {g — 2)^ and from below by LHCSewk 
constraint. The latter also provides the strongest 
upper bound on the Xilf&IXi- We also see in 
Fig. that the gluino, squark, stop and bot¬ 
tom masses are all very poorly constrained in 
our pMSSMlO analysis, though the LHCScoi con¬ 
straint forbids low masses. 

Concerning the Higgs sector, we note that the 
best-fit value for Ma lies in the multi-TeV re¬ 
gion (where its actual value is only weakly con¬ 
strained) and is therefore far in the decoupling 
region. Accordingly, the properties of the light 

®The striations in these bars reflect the n on-m onotonic 
behaviours of the function visible in Fig. |13| 


Higgs boson at about 125 GeV resemble very 
closely those of the Higgs boson of the SM. 

The first column of Table [5] lists the most im¬ 
portant contributions to the total function of 
different (groups of) constraints at the best-fit 
pMSSMlO point. The total y/ value at the best- 
fit point is y/ = 83.3, of which the largest part 
is due to the Higgs constraints evaluated using 
HiggsSignals. 

To convert the total y/ of our fit into a 
probability estimate, we calculate the contri¬ 
bution and corresponding number of degrees of 
freedom (d.o.f.) by considering only constraints 
that have signihcant contributions to our global 
y^ function in large regions of the relevant pa¬ 
rameter space. We do not include in this pro¬ 
cedure constraints from HiggsSignals, which do 
not in general vary strongly in our preferred fit 
regions (see, e.g., y/{HS) in Table Therefore, 
to calculate the y^ probability we consider in to¬ 
tal 31 constraints, which translate into 18 d.o.f 
for the pMSSMlO, 24 d.o.f. for CMSSM, 23 d.o.f. 
for NUHMl, and 22 d.o.f. NUHM2. Previous 
studies 68 showed that this definition of the y^ 
probability represents a good estimate of fit qual¬ 
ity and enables a comparison between different 
models on an equal footing. It also represents a 
reasonable approximation to the underlying ab¬ 
solute p-values of our fits. 
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Figure 11. The particle spectrum and dominant decay branching ratios at our best-fit pMSSMlO point. 
Note the near-degeneracies between XijX 2 > between the sleptons, between XsjXa emd , between 

the ql and qn, between the heavy Higgs bosons, and between the stops and bottoms, which are general 
features of our 68% CL region. On the other hand, the overall sparticle mass scales, in particular of the 
coloured sparticles, are poorly determined. 



Mfp 


Figure 12. Summary of mass ranges predicted in the pMSSMlO. The light (darker) peach shaded bars 
indicate the 95%) (68%) CL intervals, whereas the blue horizontal lines mark the values of the masses at 
the best-fit point. 
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Comparing to the ^ values for the CMSSM, 
NUHMl, and NUHM2 shown in the last three 
columns of Table [5l we see that the largest im¬ 
provement is in the contribution from (g — 2 )^, 
though there are also small improvements in 
BR(i?s —)• and the Z-pole observables. 

Overall, we see that the pMSSMlO has a prob¬ 
ability of 30.8% compared to 10.8%, 12.1% and 
11.0% for the CMSSM, NUHMl and NUHM2, 
respectively, demonstrating that the pMSSMlO 
gives a significantly better fit 

We stress, however, that these probabilities 
are only approximate and assume an underlying 
y^-distribution with no correlations between the 
observables. A more proper treatment would be 
to smear the measurements around the best-fit 
predictions, fit to these toy measurements and 
evaluate the fraction of cases in which the result¬ 
ing exceeds the observed x^ ■ We leave such an 
evaluation as a topic for future work. 


3.3. Sparticle Masses 

Fig. [^displays (from top left to bottom right) 
the one-dimensional profile likelihood functions 
for the masses of the gluino, the first- and second- 
generation squarks, the lighter stop and sbot- 
tom squarks, the lighter chargino and the lighter 
stau. In each panel the solid black line is for the 
pMSSMlO, the solid blue line for the NUHM2, 
the dashed blue line for the NUHMl and the dot¬ 
ted blue line for the CMSSM (the latter three 
lines are updated from Ref. [16| to include new 
constraints such as the LHC combined value of 
Mh 44 ). In the case of rrig, we see that signifi¬ 


cantly lower masses are allowed in the pMSSMlO 
than in the other models: > 1250 GeV at the 
68 % CL and ~ 1000 GeV at the 95% CL. We also 
see that there is a similar, though smaller, reduc¬ 
tion in the lower limit on rriq, to ~ 1500 GeV at 
the 68 % CL and ~ 1300 GeV at the 95% CL. 
The picture is more complicated for m, where 
we see structures in the one-dimensional likeli- 


^The probabilities and values of differ from those 
given in as we have updated the CMSSM and NUHM 
analyses with the most recent Higgs mass determina¬ 
tion [44] and other new information [7| |40] . We note, in 
particular, that the new BR(i?„ —> tu ^ measurement |40| 
improves the agreement with the SM and the SUSY mod¬ 
els we study. 


hood function for < 1000 GeV that reflect 
the low-mass islands in the corresponding panel 
of Fig. that are allowed at the 95% GL. In 
the bottom row of Fig. the one-dimensional 
profile likelihood functions for to-± and rrif^ in 
the pMSSM have minima at the lower mass lim¬ 
its ~ 100 GeV established at LEP, and there is 
an upper limit < 1000 GeV at the 95% CL. 
These effects are due to the (5 — 2)^ constraint 
and the choice of generation-independent slepton 
masses in the pMSSMlO. On the other hand, the 
light chargino (which is nearly degenerate in mass 
with the second lightest neutralino), has an upper 
mass limit below 500 GeV at the 95% CL. This 
would allow neutralino and chargino pair produc¬ 
tion at an 1000 GeV e+e“ collider, as we discuss 
later. 


3.4. Benchmark pMSSMlO Models 

In view of the variety of pMSSMlO parame¬ 
ters that are allowed at the 68 % CL, we con¬ 
sider in this subsection various specific bench¬ 
mark models that illustrate the range of possi¬ 
bilities. Specifically, looking at the middle panels 
of Fig. we see that a very low stop mass in 
the compressed-stop region is possible, and the 


top panels of Fig. 10 show the possibilities for a 
gluino or squark mass that is lower than at the 
best-fit point. Also, we see in the upper left panel 
of Fig. that SUSY may well appear with both 
the squark and gluino masses having lower masses 
than at the best-fit point. We investigate these 
possibilities with the benchmark points discussed 


below, whose SLHA files 37 can be downloaded 


from the MasterCode website 25 


3.4.1. Low-m^^ point 

We display in the upper left panel of Fig. M 
the spectrum at the point that minimizes 
cally within the low-m^^ (and low-mg^) 68 % CL 
region visible in the middle planes of Fig. Like 
the pMSSMlO best-fit point shown in Fig. [TT] 
this point also exhibits near-degeneracies between 
Xi,X 2 xf, between the sleptons, between 
X 3 ,X 4 and X 2 (reflected also in the fact that 
Ml ~ M 2 , as seen in the second column of 
Table|4|, and between the ql and qn. However, 
all the stops and sbottoms are light at this point. 
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Constraint d. 

o.f. 


pMSSMlO 



CMSSM 

NUHMl 

NUHM2 



best fit 

low mr 
^1 

low rrig 

low mg 

low all 
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16 


LHC 8 

1 

0.1 

1.0 

0.8 

1.0 

0.4 


- 



- 



- 


Jets-h^T 

1 

- 

- 

- 

- 

- 


2.0 


0.0 



0.5 


Mh 

1 

0.0 

0.0 

0.2 

0.2 

0.0 


0.1 



0.1 



0.4 

Mw 

1 

0.0 

0.1 

0.1 

0.5 

0.1 


0.0 


0.0 



0.4 

Bs,d g 

1 

0.2 

0.0 

0.0 

0.2 

0.0 


0.5 



0.3 



0.4 

BR (6 ^ 57 ) 

1 

0.1 

0.0 

0.2 

0.1 

0.1 


0.5 



0.0 



0.0 

BR(i3u — > TVr) 

1 

0.2 

0.2 

0.2 

0.2 

0.3 


0.2 



0.2 



0.2 


Other B physics 

5 

3.3 

3.0 

3.2 

3.3 

3.0 


3.2 



3.3 



3.3 


1 

0.1 

0.0 

0.3 

0.0 

0.1 


0.0 


0.0 



0.0 

(7p 

1 

0.0 

0.1 

0.0 

0.1 

0.5 


0.0 


0.0 



0.0 

A/H —>■ t'^t~ 

1 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 


0.0 



0.0 

Nuisance 

3 

0.0 

0.1 

0.0 

0.1 

0.8 


0.1 



0.0 



0.1 


1 

to 

1 

0.0 

0.7 

0.0 

0.0 

0.6 


9.3 

10.6 


8.4 

Z pole 

13 

16.3 

17.0 

17.1 

16.8 

16.4 

16.8 

16.5 

16.7 

Parameters 


10 -t 3 

10 -b 3 

10 + 3 

10 + 3 

10 + 3 

4 + 3 

5 + 3 

6 + 3 

xVd-o.f. 


20.5/18 

22.2/18 

22.0/18 

22.3/18 

22.2/18 

32.8/24 

31.1/23 

30.3/22 

X^ probability 


0.31 

0.22 

0.23 

0.22 

0.22 

0.11 

0.12 

0.11 

X"(HS) 

77 

62.8 

62.6 

62.8 

62.8 

62.8 

- 

- 

- 


Table 5 

Table of the total breakdowns at the pMSSMlO best-fit and low-m. 
in the CMSSM, NUHMl 


low- 


iq and low-rUg points, and 


and NUHM2 (updated from ^6 21 using in particular the current value of 
Mh^). The LHCSstop; LHCSewk and LHCScoi constraints were applied only to the pMSSMlO, whereas 
a generic jets + IpT constraint was applied to the CMSSM, NUHMl and NUHM2 ^6 21L For each set 
of constraints, the (rounded) contribution and the number of non-zero contributions is provided. The 
nuisance parameters are mt, as{Mz) and Mz- The bottom rows show the number of parameters (including 
the nuisance parameters) and the total x^/d.o.f. omitting Higgs signal rates: the latter have been calculated 


only for the pMSSMlO points, and are given separately in the last line. We also show an estimate of the 
corresponding x^ probability, which is calculated as the probability neglecting correlations between the 
observables. 


As in Fig. El the dominant decay modes are illus¬ 
trated in fifty shades of grey . The second col¬ 
umn of Table [5] lists the contributions to the total 
X^ function of different (groups of) constraints 
at this low-TO(^ pMSSMlO point. Comparing 
with the corresponding breakdown for the best- 
fit point shown in the first column of Table [H we 
see larger contributions from the LHC8 constraint 
(principally from LHCScoi) and from (g — 2)^, 
which are largely responsible for the increase in 
the total x^ to 22.2 (omitting the HiggsSignals 
contributions) and the corresponding decrease in 
the x^ probability to 0.22. However, we empha¬ 
size that this point provides a perfectly acceptable 


fit to all the constraints. 

3.4.2. Low-mq point 

We consider next a benchmark point with rel¬ 
atively low masses for the first- and second- 
generation squarks. As can be seen in the top 
right panel of Fig. the lowest value of niq that 
is allowed at the 68% CL is ~ 1500 GeV, and we 
have chosen as benchmark a point that also has 
m^o ~ 200 GeV, whose spectrum is shown in the 
upper right panel of Fig. [T^ We see there that 
the near-degeneracies between Xi, X 2 xf^ be¬ 
tween the sleptons, between X 3 ;X 4 xf, and 
between the heavy Higgs bosons are very simi- 
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Figure 13. The one-dimensional profile likelihood functions for mg, m^, m^_^, Wfi- In 

each panel the solid black line is for the pMSSMlO, the solid blue line for the NUHM2, the dashed blue 
line for the NUHMl and the dotted blue line for the CMSSM. 
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Figure 14. The particle spectra and dominant decay branching ratios at the benchmark points discussed in 
the text. Upper left panel: the low-m^^ pMSSMl 0 point, where the stops and bottoms are relatively light. 
Upper right panel: similarly for the low-my benchmark point, where all the squarks are relatively light. 
Lower left panel: similarly for the low-mg benchmark point. Lower right panel: similarly for the point 
where all squarks and the gluino masses are < 2 TeV. Note in each case the near-degeneracies between 
Xi,X 2 xt 7 between the sleptons, between X 3 JX 4 x^, between the (Jl and qn, and between the 

heavy Higgs bosons. 


lar to those at the best-fit and low-m^^ points. 
By choice, the masses of the first- and second- 
generation squarks are much lighter than at ei¬ 
ther of these points, and the third-generation 
squarks have masses intermediate between the 
best-fit and low-m^j^ points. As seen in Table [5l 
the largest part of the increase in x^ to 22.0, com¬ 
pared to the best-fit point, and the corresponding 
decrease in the x^ probability to 0.22, is again due 
to the LHC8 constraint. 


3.4.3. Low-mg point 

We consider next a benchmark point with a 
relatively low gluino mass. As can be seen in the 
top left panel of Fig.[^ our global fit requires 
my > 1250 GeV at the 68% CL. We have cho¬ 
sen as benchmark a point that has this value of 
my and also m^o ~ 200 GeV, whose spectrum is 
shown in the right panel of Fig. We see again 
the near-degeneracies within groups of MSSM 
particles, as for the benchmark points considered 
previously. We see a clear hierarchy of masses 
between the groups of strongly-interacting spar- 
ticles, with the third-generation sparticles being 
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much heavier than those of the first and second 
generation, which are in turn much heavier than 
the gluino. Again as seen in Table O the largest 
part of the increase in ~x^ —)■ 22.3 compared to 
the best-fit point is again due to the LHCScoi con¬ 
straint, with increases also from LHCSewk, and 
Mw- The total ^ probability of 21.7% is com¬ 
parable to those of the low-m^^ and -m^ points. 

3.4.4. Point with squark and gluino masses 
below 2 TeV 

Finally, we display in the lower right panel 
of Fig. the spectrum at a point from near 
the turning-point in Fig. which can be re¬ 
garded as a ‘compromise’ between the two pre¬ 
vious benchmarks where the gluino and all the 
squarks (including those in the third generation) 
have masses < 2 TeV, as do the heavy AjH Higgs 
bosons. Like the previous pMSSMlO benchmark 
points, this point also exhibits near-degeneracies 
between Xi, X 2 xf, between the sleptons, be¬ 
tween X31X4 xt j s-nd between the and 
qR. In addition, the sbottom squarks are also 
nearly degenerate, whereas the stops exhibit a 
greater mass splitting, due to the mt-dependence 
in the off-diagonal stop mass matrix elements. 
The contributions to the total x^ function of dif¬ 
ferent (groups of) constraints at this low-mass 
pMSSMlO point are shown in the fifth column 
of Table 1 Comparing with the corresponding 
breakdown for the best-fit point, we see larger 
contributions from LHCScoi, and (g — 2)^. All 
these contributions to x^ are < 1, but they do 
suggest possibilities for SUSY discovery in jets + 
.^Tsearches early during LHC Run 2, and in up¬ 
coming direct dark matter detection experiments. 
The total x^ = 22.2 at this point, and the corre¬ 
sponding x^ probability is 22.4%. 


3.5. The Anomalous Magnetic Dipole Mo¬ 
ment of the Muon (g — 2)^ 

It is well-known that there is a discrepancy of 


3.5cr between the measured value of (5—2)^ 15 


and the value predicted in the SM [70[|71| . Size¬ 
able contributions to (g — 2)^ from SUSY can 
occur when smuons, charginos and the lightest 
neutralino have masses of 0(100 GeV). It is 
known from previous analyses of the CMSSM, 


NUHMl and NUHM2 16 21 that in these mod¬ 
els there is tension between SUSY interpretations 
of the discrepancy in the anomalous magnetic 
dipole moment of the muon (g — 2)p (which favour 
lower electroweak sparticle masses) and LHC con¬ 
straints from direct searches for sparticles and 
the measured value of the lightest Higgs boson 
(which favour higher coloured sparticle masses). 
This tension arises from the universality relations 
imposed in these models at the GUT scale be¬ 
tween the soft SUSY-breaking contributions to 
the masses of the strongly- and electroweakly- 
interacting sparticles. In the pMSSMlO there are 
no such assumptions, and thus one might hope to 
resolve this tension. 

This point is apparent in Fig. |15[ where we dis¬ 
play in the left panel the contributions from 
(g — 2)p to the global x^ functions of our fits to 


the CMSSM (blue dotted line), the NUHMl (blue 
dashed line), the NUHM2 (blue solid line) and the 
pMSSMlO (black solid line), as well as the exper¬ 
imental likelihood function that we assume (solid 
red line). We see that the pMSSMlO is able to 
fit (g — 2)^ perfectly with ~ 0 at the best- 
fit point, whereas the other “universal” models 
exhibit contributions Ay^ > 9 from (g — 2)^. 

We display in the right panel of Fig. [T^ 
the impact on the global y^ as a function of 
(g — 2)^ of implementing the LHC constraints 
on electroweakly-interacting sparticles using the 
LHCSewk method described earlier (which, as 
we have shown, provides a reasonably accurate 
as well as computationally economical represen¬ 
tation of the LHC8 constraints on electroweakly- 
interacting sparticles). The solid line is the 
global x^ function with the LHCSewk constraint 
included, and the dashed line when they are 
omitted. The minimum value of y^ increases 
from 82.6 to 83.3, and the value of (g — 2)^ 
at the minimum is essentially unchanged. We 
conclude that the impacts of the LHC searches 
for electroweakly-interacting particles are limited, 
and the pMSSMlO resolution of the (g — 2)^ puz¬ 
zle survives the LHC electroweak constraints with 
flying colours. 

The left panel of Fig. displays the two- 
dimensional profile likelihood function in the 
(mpjj,m^o) plane, with the solid (dashed) 
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Figure 15. Profile likelihoods for the SUSY contribution to {g — 2)^. The left panel shows the Ax^ 
contributions from (g — 2)^ to the global likelihood functions of our fits to the CMSSM (blue dotted line), 
the NUHMl (blue dashed line), the NUHM2 (blue solid line) and the pMSSMlO (black solid line), as well 
as the experimental likelihood function that we assume (solid red line). The right panel displays the global 
function calculated without (dashed line) and with (solid line) the contribution of the electroweakly- 
interacting sparticle searches implemented via LHCSewk- 


red/blue contours denoting the Ax^ = 2.30/5.99 
level contours for the case where we do (not) 
apply the LHC 8 constraints, respectively, and 
the green filled and empty stars indicating the 
corresponding best-fit points Qualitatively, 
this plane is quite similar to the corresponding 
(m?j, m^o) plane shown in the bottom right panel 
of Fig. [l^ though we note, e.g., that the best-fit 
value of is ^ 100 GeV larger than the best- 
fit value of mfj. This feature is apparent also 
when one compares the right panel of Fig. |16[ 
which displays the one-dimensional profile like¬ 
lihood function for with the corresponding 
plot for rUf^ in the bottom right panel of Fig. [TSl 
In both cases, the one-dimensional profile like¬ 
lihood function in the pMSSMlO is shown as a 
solid black line, that in the NUHM2 as a solid 
blue line, that in the NUHMl as a dashed blue 
line, and that in the CMSSM as a dotted blue 
line. 


®We do not show the corresponding results for the UL, 
which are very similar. 


3.6. Interplay of the LHCSewk? (g — 2)^ 
and Dark Matter Constraints 

The 68 % and 95% CL regions in the 
(m^^,tan/3) plane before (dashed lines) and af¬ 
ter (solid lines) implementation of the LHC 8 and 
other constraints are displayed in Fig. We 
see that the lowest values of tan/3 receive a x^ 
penalty, which is due to a combination of different 
effects. In particular, the LHC8ewk constraint 
disfavours lower values of ^ which, in com¬ 
bination with {g — 2)^, results in a x^ penalty 
for tan /3 < 10. Because we impose slepton 
mass universality in the pMSSMlO, stau masses 
are also pushed to higher values. In this way 
the LHC8ewk constraints eliminate pMSSMlO 
models with a Bino-like LSP and small for 
which stau coannihilation and t-channel slepton 
exchanges brought the relic LSP density into 
the allowed range. The remaining models with 
tan ,3 < 30 then fall foul of the LUX upper limit 
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on (Tp^, because the LSP has a substantial 
Higgsino component, which enhances Upb 


The 


overall combined effect of the LHC 8 ewk, (g — 2 )^ 
and dark matter constraints is to prefer values of 
tan/3 between about 15 and 45 at the 68 % CL, 
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★ - - pMSSMlO w LHC8: best fit, Icr, 2 ct 


•sS- - - pMSSMlO w/o LHC8: best fit, Icr, 2cr 




Figure 16. Left panel: The two-dimensional profile likelihood function for versus the lightest neu- 
tralino mass m^^o. The solid (dashed) red/blue contours denote the Ax^ = 2.30/5.99 level contours for 
the case where we do (not) apply the LHC8 constraints, respectively, and the green filled and empty stars 
indicate the corresponding best-fit points. Right panel: The one-dimensional profile likelihood function for 
in the pMSSMlO (solid black line), the NUHM2, (solid blue line), the NUHMl (dashed blue line) 
and the CMSSM (dotted blue line). 


though tan/3 values below 10 are still allowed at 
the 95% CL. We note that this feature is an effect 
of the choice of a single slepton mass scale, which 
could be avoided in more general versions of the 
pMSSM. 

3.7. Higgs Physics 

Fig. displays one-dimensional profile likeli¬ 
hood for Mh when the LHC constraints are ap¬ 
plied. We see that the likelihood for Mh in the 
pMSSMlO (black line) is very similar to the ex¬ 
perimental value smeared by the theoretical un¬ 
certainty in the FeynHiggs calculation of Mh for 
specific values of the MSSM input parameters. 

The left panel of Fig. displays the two- 
dimensional profile likelihood function in the 
{Ma, tan fi) plane. As before, the solid (dashed) 
red/blue contours denote the = 2.30/5.99 
level contours for the case where we do (not) 
apply the LHC8 constraints, respectively, and 
the green filled and empty stars indicate the 


corresponding best-fit points. Comparing the 
dashed and solid 68% contours, we see that 
lower values of tan/3 are disfavoured at the 68% 
CL by the combination of LHC8 ewK) {g — 2)^ 
and Dark Matter constraints, as discussed in 
the previous subsection. Those constraints, in 
combination with the choice of a single slep¬ 
ton mass scale for all three generations, lead 
to limits of Ma ^ 1000(500) GeV at the 68 
(95)% CL, whereas otherwise low CP-odd Higgs 
boson masses down to Ma ~ 500(350) GeV 
would be found in the 68 (95)% CL area. 

The right panel of Fig. displays the corre¬ 
sponding one-dimensional profile likelihood func¬ 
tion for Ma- as before, the solid black line is for 
the pMSSMlO, the solid blue line for the NUHM2, 
the dashed blue line for the NUHMl and the dot¬ 
ted blue line for the CMSSM. Lower Ma values 
for tan /3 < 30 are in particular disfavoured by the 
LUX and other limits, as discussed in the previ¬ 
ous subsection. 
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★ - - pMSSMlO w LHC8: best fit, la, 2a 


* - - pMSSMlO w/o LHC8: best fit, la, 2a 



Figure 17. The 68% and 95% CL regions in the (m^^,tan/3) plane before (dashed lines) and after (solid 
lines) implementation of the LHC8 and other constraints. 



Figure 18. The one-dimensional profile likelihood function for M^: the solid black line is for the 
pMSSMlO, the solid blue line for the NUHM2, the dashed blue line for the NUHMl, the dotted blue 
line for the CMSSM, and the red line is the penalty from the experimental measurement of Mh with 
the assumed theoretical uncertainty o/1.5 GeV. 
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★ - - pMSSMlO w LHC8: best fit, Icr, 2 ct 

* - - pMSSMlO w/o LHC8; best fit, Icr, 2cr 




Figure 19. Left panel: The two-dimensional profile likelihood function for Ma versus tan/3. The solid 
(dashed) red/blue contours denote the Ax^ = 2.30/5.99 level contours for the case where we do (not) 
apply the LHC constraints, respectively, and the green filled and empty stars indicate the corresponding 
best-fit points. Right panel: The one-dimensional profile likelihood function for Ma-' the solid black line 
is for the pMSSMlO, the solid blue line for the NUHM2, the dashed blue line for the NUHMl and the 
dotted blue line for the CMSSM. 


3.8. BR(Bs —>■ pi~) Decay 

We display as a black line in Fig. the pro¬ 
file likelihood in the pMSSMlO for the ratio of 
BR(i?s —) to the SM value. This can be 
compared with the penalty from the experi¬ 
mental constraint on BR(i3s —>■ pfipT), which is 
shown as a red line. It is interesting to note that 
in the pMSSMlO both enhancement and suppres¬ 
sion are possible, as opposed to the CMSSM, the 
NUHMl and the NUHM2 16pT , in which a sup¬ 
pression was not possible and only an enhance¬ 
ment was allowed. This comes about because the 
extra parameters in the pMSSMlO make possible 
some negative interference between the SM and 
SUSY amplitudes, which is not possible in the 
other models when the various other constraints 
are implemented. 


3.9. Direct Dark Matter Detection 

The left panel of Fig. displays the one¬ 
dimensional profile likelihood in the pMSSMlO 
for rn^o with the same colour coding as in Fig. |13| 


We see that, in contrast to the other models, the 
pMSSMlO favours a low mass for the xf), driven 
again by the (g — 2)^ constraint. The right panel 
of Fig. displays the two-dimensional profile 
likelihood for the lightest neutralino mass versus 
the spin-independent cross-section, where the red 
and blue contours show the 68% and 95% CL lev¬ 
els respectively. The region that is excluded by 
LUX 42 and XENONIOO 43 is shaded green, 
whereas the ‘floor’ below which the background 
from atmospheric neutrinos dominates is shaded 
yellow 73 . The low-mass vertical 95% CL strips 
are due to points where the relic LSP density is 
brought into the cosmological range by annihila¬ 
tions through direct-channel Z and h poles. 

It is interesting to note that the pMSSMIO fit 
prefers rather high values of the spin-independent 
cross section after application of the LHC8 con¬ 
straints: lower values could be reached for a 
Bino-like LSP, but the dark matter density con¬ 
straint would then require stau coannihilation 
and t-channel slepton exchange, which are, how- 
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Figure 20. The one-dimensional profile likelihood functions for —>■ pf'pL~) relative to the SM value. 

The solid black line is for the pMSSMlO, the solid blue line for the NUHM2, the dashed blue line for the 
NUHMl, the dotted blue line for the CMSSM, and the red line is the penalty from the experimental 
constraint. 


★ - - pMSSMlO w LHC8: best fit, la, 2a 

* - - pMSSMlO w/o LHC8: best fit, IcT, 2(7 




Figure 21. Left panel: The one-dimensional profile likelihood in the pMSSMlO for (black line), 
compared with the NUHM2, the NUHMl and the CMSSM (solid, dashed and dotted blue lines, respec¬ 
tively) .Right panel: The two-dimensional profile likelihood function in the pMSSM in the {m.^o,ap^)-plane, 
.showing the regions excluded by the XENONIOO and LUX experiments (shaded green), the neutrino ‘floor’ 
(shaded yellow), and the prospective sensitivity of the LZ experiment (purple) f7^y. 


ever, disfavoured by the combined effects of the 
LHCSewk and {g — 2)^ constraints. Our best-fit 


region is close to the present experimental up¬ 
per limit on (7®^ 42 , and consequently within 
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reach of future direct detection experiments such 
as LZ 72 , as indicated by the magenta line in 
the right panel of Fig. On the other hand, 
we note that before applying these constraints, 
and even afterwards at the 95% CL, there are 
values for that go far below this neutrino 
‘floor’, highlighting the complementarity of direct 
detection experiments and searches at the LHC. 
Since these very low values of are due to can¬ 
cellations between different contributions to the 
spin-independent scattering matrix element, one 
may ask whether the spin-independent cross sec¬ 
tions on proton and neutron targets could be very 
different when this cancellation occurs. More 
specifically, one may wonder whether, for mod¬ 
els in which is below the neutrino ‘floor’, the 
cross-section for scattering on a neutron tar¬ 
get may be less suppressed, perhaps remaining 
above the neutrino ‘floor’ ? As we see in Fig. 
the spin-independent cross sections on proton and 
neutron targets are generally very similar when 
> 10 “"*^ cm^, but may indeed be quite dif¬ 
ferent when (Tp^ < 10 “^® cm®, which is approx¬ 
imately the lowest level of the neutrino ‘floor’, 
whose height varies as seen in the right panel 
of Fig. Points coloured black (green) [blue] 
{red} have both and above the neutrino 
‘floor’ shown in the right panel of Fig. ( 17 ®^ 
below and above) [tr®^ above and cr®^ below] 
{ct®^ and CT®^ both below}. We see that there 
is a significant population of models whose spin- 
independent scattering cross sections on protons 
and neutrons are both below the ‘floor’ (indicated 
in red), so there is no ‘no-lose’ theorem for dark 
matter scattering in the pMSSMlO 


4. Extrapolation to High Scales 

In our analysis of the pMSSMlO we have not 
imposed any restriction on the possible extrap¬ 
olation of the (purely phenomenological) soft 

®We do not include in the right panel of Fig. |21| and in 
Fig. the contributions of loop-induced scattering off 
gluons [74] . In general, these contributions are relatively 
small [75| , but they would also shift slightly the param¬ 
eters of the models exhibiting strong cancellations in cTp^ 
and cr^^. We thank N. Nagata for discussions on these 
points. 


SUSY-breaking parameters to high scales using 
the renormalization-group equations (RGBs). In 
many cases, one could expect that renormali¬ 
sation by the gaugino masses may drive some 
soft supersymmetry-breaking sfermion masses- 
squared TOq to negative values at high-energy 
scales 76 . This raises cosmological issues that 
have been studied, for example, in 77 , and such 
scenarios do not necessarily lead to an unaccept¬ 
able evolution of the Universe. However, it is in¬ 
teresting to study the implications of requiring 
TOq > 0. We emphasise that this cut reduces the 
data set significantly, and one may anticipate that 
part of the parameter space would be recovered 
in a dedicated scan. Nevertheless, we expect that 
the main features discussed here would be present 
also in a more complete scan. 

Fig. [^displays the two-dimensional likelihood 
functions in some relevant sparticle mass planes. 
In each panel, the red (blue) lines are the 68 % 
(95%) CL contours, the solid (dashed) lines be¬ 
ing after (before) a cut requiring TOq > 0 for 
all the sleptons and squarks at the GUT scale 
^ 2 • 10^® GeV. The upper left panel shows 
the (TOq,TOg) plane, and can be compared with 
Fig-i (upper left plot). We see that the pri¬ 
mary impact of the anti-tachyon cut is to remove 
all models above a diagonal line where the neg¬ 
ative renormalisation by M 3 drives the squark 
masses-squared negative at the GUT scale. The 
upper right panel of Fig. shows the impact 
of the anti-tachyon cut on the (TO(^,to^o) plane, 
which can be compared with the middle left 
panel of Fig. Here the most obvious impact 
is to remove the compressed stop region where 
TO(^ —m^o < TOtF% The lower left panel of Fig. 
shows the impacFin the (to-±, Tn^o) plane, where 
we see that many models with small m~± — m^o 

Xl A1 

survive the anti-tachyon cut. The anti-tachyon 

cut leads to a more pronounced preference for 
small values of m^o and in particular m-±. Fi- 
nally, the lower right panel of Fig. [^displays the 
(to^^ , rrij^o) plane, where we see that the anti- 
tachyon cut has very little effect, except to re¬ 
move some points with small to^^ — In 


^‘^The 68% CL region extends to much larger values of 
where larger values of are also found. 
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proton log((j®Vcm^) 


Figure 22. Scatter plot of the cross sections for spin-independent scattering on a proton target (horizontal 
axis) and on a neutron target (vertical axis) obtained from a sampling of pMSSMlO points within the 
95% CL region. The diagonal dashed line corresponds to equal spin-independent cross sections on proton 
and neutron targets. The colour-coding distinguishes between points with either and/or cr®^ above or 
below the neutino ‘floor’ seen in\21\ 


particular, the best-fit values of and rn^^o 

are little changed, and the mitigation of the 
{g — 2)^ anomaly in the pMSSMlO survives the 
anti-tachyon cut. However, we repeat that this 
cut may even not be necessary 

Fig. [24| shows the impacts of the optional anti- 
tachyon cut on the one-dimensional profile likeli¬ 
hood functions for mg, mz , m^o, m~± and 

(from top left to bottom right). We see that 
the function for the gluino mass is little af¬ 
fected, whereas points with low mq are system¬ 
atically removed, as one might expect from en¬ 
forcing TOq > 0. These effects can also be seen in 
the upper left panel of Fig. As one would ex¬ 
pect from the upper right panel of Fig. |23[ points 
with low are also removed by the anti-tachyon 
cut, and the best-fit value of mi^ is increased by 
~ 1 TeV. As seen in the middle right panel of 
Fig. the one-dimensional likelihood function 
for r7T,<;o is little affected, whereas that for m-± 

AT ' Yn 


is squeezed strongly. These effects reflect the 
behaviour in the (m~± , m^o ) plane seen in the 
lower left panel of Fig. where the favoured 
points lie in a narrow // — Xi coannihilation 
strip. These points have Mi ~ M 2 at the elec- 
troweak scale, leading to the potential observabil¬ 
ity of neutralino/chargino pair production at an 
e“'"e“ collider with a centre-of-mass energy be¬ 
low 1000 GeV, as we discuss later. Finally, we 
see in the bottom right panel of Fig. that the 
likelihood function for is little affected by 

the anti-tachyon cut, apart from the removal of 
some low-mass points as seen already in the lower 
right panel of Fig. However, as already com¬ 
mented, the removal of these points does not pre¬ 
vent the pMSSMlO from addressing successfully 
the {g — 2)^ problem. 

As a final topic in this Section, we dis¬ 
cuss the departures from universality of the soft 
supersymmetry-breaking parameters in the sam- 
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Figure 23. The impacts of the optional anti-tachyon cut on the two-dimensional profile likelihood functions 
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and ) planes. In each panel the solid (dashed) red/blue 


contours denote the Ax^ = 2.30/5.99 level contours for the case where we do (not) apply the anti-tachyon 
constraint, respectively. The green filled and empty stars indicate the corresponding best-fit points. 


pie that would survive the anti-tachyon cut. 
Fig. shows a plane of the root-mean-squared 
deviations from gaugino- and sfermion-mass uni¬ 
versality, defined by 




N 





( 4 ) 


where the denote, respectively, the various 
gaugino mass parameters and the square roots 
of the (positive) squark and slepton mg pa¬ 
rameters in the pMSSMlO at the GUT scale, 
and rh denotes their respective averages. Exact 
unification of the gaugino (sfermion) masses is 
achieved when ctm (fm) vanishes. We see that 
sfermion-mass universality is quite strongly vio¬ 


lated, and gaugino-mass universality is also dis¬ 
favoured, though still possible at the 95% CL. As 
we have already commented, the favoured points 
in the narrow — x^ coannihilation strip must 
have near-degenerate Xi and X 2 and hence M 2 — 
Ml at the SUSY-breaking scale, corresponding 
to a breakdown of universality by a factor ~ 2 at 
the GUT scale, i.e. MiI^Mquji) ^ 2 Af 2 (Af(juT)- 
As can also be inferred by comparing the top left 
and middle right panels of Fig. a violation 
of GUT-scale M 3 — Mi universality is also sug¬ 
gested. Thus, refined future fits based on more 
data might lead to a preference for some different 
scenario for unification. 
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Figure 24. The impacts of the optional anti-tachyon cut on the one-dimensional profile likelihood functions 
for mg, mq, mf_^, rn^o, iti^± and In each panel the solid (dashed) lines are for the cases where we 

do (not) apply the anti-tachyon constraint, respectively. 
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Figure 25. Two-dimensional likelihood function in the plane of the root-mean-square deviations from 
sfermion- and gaugino-mass universality, cFm O'^d aM, defined in the text. 


5. Prospects for Sparticle Detection in Fu¬ 
ture LHC Runs 

At the time of writing, the LHC is starting 
Run 2, taking data at 13 TeV, and it is expected 
that an integrated luminosity of 300 fb~^ will be 
collected by the early 2020s. There are also plans 
for a subsequent high-luminosity upgrade to ac¬ 
cumulate 3000 fb“^. In this Section we describe 
some prospects for future direct LHC searches for 
sparticles by ATLAS and CMS that follow from 
our analysis of the pMSSMlO. 

With the increase of the LHC centre-of-mass 
energy from 8 TeV to 13 TeV for Run 2, there 
will be large increases in the reaches for high- 
mass sparticle states. As shown in Fig[T^ gluino 
masses ~ 1.25 TeV (top left panel) and first- 
and second-generation squark masses ~ 1.5 TeV 
(top right panel) are within our 68% CL region. 
These masses will be probed by ATLAS and CMS 
with just a few fb~^ of data, demonstrating that 
already in an early phase of Run 2 the discov¬ 
ery of SUSY might well be possible. For third- 
generation squarks, it is important to point out 


that besides masses of ^ 800 GeV for U (middle 
left panel) and ^ 1 TeV for sbottoms (middle left 
panel), we also find in our 95% CL region masses 
that are ^ 200 to 600 GeV in the compressed stop 
region and ^ 500 GeV for sbottoms. These re¬ 
gions have not been excluded by the LHC searches 
so far, but should become partly accessible in 
the first years of 13-TeV operation. As we com¬ 
ment later, in the cases of compressed-spectrum 
charginos (bottom left panel) and sleptons (bot¬ 
tom right panel) comprehensive coverage of the 
preferred parameter space in the pMSSMlO by 
the LHC experiments will be challenging. How¬ 
ever, depending on the decay modes of the elec- 
troweakly produced sparticles, early discovery at 
13 TeV might also be possible. 

Turning to the long-term prospects for the 
LHC, the ATLAS Collaboration has made physics 
studies that explore the discovery and exclu¬ 
sion reach of ATLAS with 300 and 3000 fb“^ 
at 14 TeV: see Fig. 13 of [^. In Fig. [M| we 
display in the {mq,mg) plane our 68% (95%) 
CL contours in red (blue) as well as the esti- 



37 


* - - pMSSMlO W LHC8: best fit, la-, 2a 



Figure 26. The (mg,mg) plane with our 68 and 95% CL contours shown as solid red and blue lines, 
respectively, and the best-fit point as a green star. Also shown as solid (dashed) magenta lines are the 
estimated ATLAS sensitivities for 5-a discovery (95% CL exclusion) of SUSY via the generidfr search 
with 300 fb~^ at 14 TeV. 


mated 5-a discovery (95% CL^ exclusion) sensi¬ 
tivity with 300 fb“^ as solid (dashed) magenta 
contours This shows that a substantial region 
of our preferred parameter space, including our 
best-fit point, is within reach of future LHC runs. 
However, we recall that the position of our best- 
fit point in the (rnq,mfi) plane is rather poorly 
determined. 

In the following we revisit the mass planes of 
Fig. assessing carefully the decay modes of the 
respective SUSY particles. A recurring theme is 
that the and y® nearly degenerate in mass 
with in the 68% CL region, so that squarks 
and sleptons decay via or X 2 large fractions 
of the preferred parameter space. This general 
scenario is consistently indicated using pale blue 
shading. 

With this in mind we turn to Fig.[^ where 
in the upper left panel we explore the possi- 

^^The 5-a discovery contour for 3000 fb“^ is almost coin¬ 
cident with the 95% exclusion contour for 300 fb“^. 


ble future LHC sensitivity to direct stop pro¬ 
duction in the compressed-spectrum region. As 
previously, our present 68% (95%) CL contours 
are shown in red (blue). The colour shadings 
code the regions where the corresponding branch¬ 
ing ratio, shown in the legend, exceeds 50% for 
the point at each location that minimises the 
X^ function over the remaining parameters, and 
the thin diagonal dashed black lines correspond 
to Am = mf^ — m^o = 0,Mw -\- mb and m^. 
The solid dashed black lines show the projected 
LHC 95% CLg exclusion sensitivities for U —f 
decays with 300 fb~ 


79 


(similar sensitivity is 
found in this region with 3000 fb“^). These do 
not cover the case of a compressed spectrum re¬ 
gion, which includes the 95% CL region where the 
dominant ii decays are to x^b. Here we rescale 
from the present 95% CLg limit from the dibot¬ 
tom analysis, assuming that ~i^xi ^ GeV 
and using the Collider Reach tool 80 to rescale 
the production cross-section, and assume that fu- 
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Figure 27. Upper left panel: The plane with our 68 and 95% CL contours shown as solid 

red and blue lines, respectively, as well as coloured regions where the indicated branching ratios exceed 
50%. The projected LHC sensitivity with 300 fb“^ for ii —f Xi + ^ decays is shown as a thick black 
line, and the corresponding sensitivity for ii —f xfb decays is shown as a pale blue dashed line. Upper 
right panel: The {m.^±,m.^o) plane with our 68 and 95% CL contours shown as solid red and blue lines, 
respectively. The shadings indicate where the branching ratios exceed 50%. Also shown as solid (dashed) 
yellow/orange/purple lines are the projected LHC 95% CLg exclusion reaches for associated xf and X 2 
production with decays via WjZ/WIh/lLjvi^/Ti^jvT^ with 300 (3000) fb~^ of data if these decays are 
dominant. Lower left panel: The plane with our 68 and 95% CL contours shown as solid red 

and blue lines, respectively, with pale blue shading showing also where the branching ratio for fiR —>■ /xx? 
is dominant, typically > 90%. The solid (dashed) pale blue lines show our estimates of the LHC 95% 
exclusion reach with 300 (3000) fb“^. Lower right panel: Similarly for the plane, displaying 

the regions where the fiL —f MXiJ or decay modes have branching ratios exceeding 

50%. The red lines indicate the 95% exclusion reach with 300 (3000) fb“^ if Ol —f UXi were dominant, 
but are also indicative for the decay into xf /X 2 ii^ot have masses nearly degenerate with Xi. 
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Figure 28. As in Fig. 21 but for g decays (left panel) and q decays (right panel). The pale blue solid 
(dashed) lines show the estimated sensitivities with 300 fh~^ (3000 ) for (left panel) g 

and (right panel) pl ^ 9 X 2 
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ture LHC searches maintain the same search per¬ 
formance, i.e., the same signal yield after the 
event selection as present searches. We see that 
a search with 300 fb“^ of data (pale blue line) 
would already cover part of the 95% CL region 
in the compressed-spectrum region, and the esti¬ 
mate for 3000 fb~^ is similar. 

In the upper right panel of Fig.[^ we ex¬ 
plore the possible future LHC sensitivity in the 
(m~±,TOj^o) plane where our current 68% 
(95%) CL contours are again shown in red (blue), 
and the best-fit point is indicated by a green star, 
and the thin diagonal dashed black lines corre¬ 
spond to Am = — m^o = 0,Mz and Mh. 

We use colour-coding to display points in the 
(m~±,TOj^o) plane where the following xf/X 2 de¬ 
cay modes have branching ratios > 50%: via vir¬ 
tual bosons ff'Xi/xl ffXi (pale blue), 

via on-shell bosons xf bFXi/X 2 ^Xi 
low) or xf —t IFXi/X 2 ^Xi (orange), via slep- 


recall that in the part of this region favoured at the 
68% CL the x? ~ X 2 X 2 ~ xf mass differences are 
both small. 


tons Xi where = 

e,pL) (red) and xf VrfL{TVr)/X 2 TfLi^r^r) 
(purple), whereas points with no branching ratios 
exceeding 50% are coloured grey. The ATLAS 
Collaboration has made available projections of 
its sensitivities for some relevant searches for asso¬ 
ciated xf and X 2 production with 300 (3000) fb“^ 
of data 81 , which are also shown in the upper 
right plane of Fig. as solid (dashed) contours 
in the same colours as the relevant decay modes: 
yellow for xfx 2 via WZ, orange for xfx® via 
Wh, red for xf X 2 via where {£ = e,p), and 
purple for xf X 2 via t\Ivt- With 300 fb“^ of data 
the Wh search should already cover essentially all 
of the 95% CL island with m^o < 80 GeV, where 
these branching ratios exceed 50% However, 
even with 3000 fb~^ of data these searches would 
have limited impact on the other 95% CL regions, 
since there the branching ratios for these decays 
are typically small. More importantly, they would 


^^However, it should be kept in mind that the projected 
exclusion regions always assume the relevant branching 
ratios to be equal to one, and real exclusion bounds with 
“mixed” branching ratios would in general be different. 
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have no impact on the 68% CL region. The 
most relevant searches in the region with near¬ 
degenerate xf , X 2 and X? would be in hadronic 
final states sensitive to Xi —t ff'Xi/X 2 ffXi 
and xf —>■ VTTL{ji/r)/x 2 Searches 

for compressed charginos/neutralinos have been 
explored in [82| , where some sensitivity was found 
up to m~± /m^o < 300 GeV (after 3000 fb“^), al¬ 
though for very small mass differences (to-± ~ 
TU-^o) — rn^o < 20 GeV and with optimistic as¬ 
sumptions about the possible systematic uncer¬ 
tainties. 

The lower left panel of Fig. provides infor¬ 
mation about the dominant branching ratio for 
the /xij in the favoured region of the , m^o ) 
plane: as in the previous panels, our current 
68% (95%) GL contours are in red (blue). We 
see that the branching ratio for —>■ /iXi ex¬ 

ceeds 50% in all of the 95% CL region. We also 
show projections of the possible future sensitiv¬ 
ities of the LHC with 300 (3000) fb“^ of data 
to jxR —)■ /rxi decay as solid (dashed) pale blue 
lines. These projections were obtained via the 
following steps: 1) the present LHC 95% CL^ 
limit for large J m^o was rescaled using the 
Collider Reach tool [80| to estimate the flu pro¬ 
duction cross-section, and assuming that future 
LHC searches maintain the same search perfor¬ 
mance as present searches, and 2) we assumed 
that the shapes of the future sensitivity curves for 
other values of mpj^/mp^o would be the same as for 
the current searches. We see that with 300 fb“^ 
the LHC would already explore a substantial part 
of the current 68% CL region in the , m^o) 
plane, and that most of the 95% CL region could 
be explored with 3000 fb~^ but missing a narrow 
band where is small. 

The favoured region of the {mp^,mp.a) plane, 
shown in the lower right panel of Fig. looks 
similar, but the dominant decay modes are more 
varied: any of the decays /11 —>■ ^J-Xl, /t^X 2 

or Vfj,xt/HXi have branching ratios exceeding 
50%. However, the fiL ^ MXi decay mode dom¬ 
inates only when is very small. We 

have used the same approach as used above for 
projecting the fin sensitivity to the ptXi decay 
mode also to estimate the future jiL sensitivity. 


as shown by the solid (300 fb“^) and dashed 
(3000 fb“^) deep red lines. Although this pro¬ 
jection assumes decays directly into Xij it may 
have similar sensitivity to the decay into Xi^/X 2 - 
displays the corresponding (TOg,TO^o) 


Fig. 
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and (ruq,771^0) planes. In the g case (left 
panel), we see that a number of different decay 
modes may have a branching ratio exceeding 50%: 
via off-shell first- and second-generation squarks 
into charginos/neutralinos g —>■ qq'xf/qqX 2 
(pale blue), through off-shell third-generation 
squarks into (heavier) charginos/neutralinos g —>■ 
3 4 (green), into first- and second- 
generation squarks g ^ qq (yellow) or third- 
generation squarks g —)■ qsqs (orange). How¬ 
ever, none of these decay modes may have a 
branching ratio exceeding 50% (grey). In the q 
case (right panel), the decays of q ^ q'xt/qX 2 
(pale blue) are usually dominant, particularly at 
lower masses, though in some cases decays into 
gluinos qL ^ qg and at low m^o decays into heav¬ 
ier neutralinos/charginos qL —t qXi/q'xt (pale 
green) are dominant. The ATLAS Collabora¬ 
tion has presented projected exclusion limits for 
9 qqxi and q —)■ gxi simplified models with 
300 (3000) fb-i of data at 14 TeV in [sT]. Re¬ 


calling that rn^± (m^o) ~ m^o in our 68% CL 
region, these simplified model limits may be ap¬ 


plicable in the pale blue regions in Fig. 28 where 
9 qq'Xi/qqX 2 iq q'Xi/9X2) dominate the 
gluino (squark) decay modes. We overlay these 
projected limits in Fig. 28 The solid (dashed) 


curves corresponds to the 300 (3000) fb“^ data. 
In the {rrig, plane we can see that a large 

part of our 68% and 95% CL regions can be 
probed with 300 (3000) fb“^ data. Indeed, our 
best-fit point lies on the projected limit for 3000 
fb“^. We also see in the (m,, m^o) plane that 
some parts of our 68% and 95% CL regions can 
be explored with 300 (3000) fb“^ data, although 


the projected limit presented in 81 with 3000 
fb“^ data does not reach our best-fit point. 

Finally, we turn the prospects for discov¬ 
ery with 300 fb”^ of our benchmark points, 
starting with our global best-fit point (Fig. 0, 
which is just inside the reach of generic ]^t 
searches (see Fig. [26|), well in reach for the 
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slepton searches (lower panels of Fig. [^, and 
even potentially within reach of the compressed- 


chargino/neutralino searches, as discussed in 82 


due to its small mass splitting rn^± — m^o ~ 


20 GeV. 

As for our local best-fit point in the low-m^^ 
region (see the upper left panel of Fig.[^, it 
lies just within the reach of future searches in 
the compressed-stop region (upper left panel of 
Fig. [27|), a s well as slepton searches (lower panels 
of Fig. |27D, but would be difficult to access via 
chargino/neutralino searches, because of the low 
mass splittings seen in Fig. [T4| and the relatively 
high m^o ~ 300 GeV. The relatively large con¬ 
tribution of the LHGScoi constraint for this point, 
seen in the third row of Table [5l indicates that 
this point may be accessible via jets + X + I^t 
searches early in Run 2. 

In the cases of the low-m^ and/or -rrig points, 
by construction these points could also be dis¬ 
covered early in Run 2 of the LHC, since they lie 
very close to the current 68 % CL boundary in the 
{mq,mg) plane shown in Fig. |26| This feature is 
also indicated by the significant contributions to 
the global functions for these points that can 
also be seen the third row of Table [H 


6. Prospects for Sparticle Detection at a 
Future e+e“ Collider 

Fig. displays the one-dimensional x^ func¬ 
tions for the lowest particle pair- and associ¬ 
ated chargino and neutralino production thresh¬ 
olds in e+e” annihilation in the pMSSMlO 
(black), compared with their counterparts in the 
CMSSM (dotted blue), NUHMl (dashed blue) 
and NUHM2 (solid blue). In the cases of XiXi 
(upper left panel), X 1 X 2 (upper right panel) and 
xf xf (lower right panel) production, we see that 
the minima of the x^ functions in the pMSSMlO 
lie within reach of an e“'"e“ collider with centre- 
of-mass energy 500 GeV, and that threshold loca¬ 
tions favoured by < 3 would be within reach 
of a 1000 GeV collider, whereas no upper limit 
can be established at the 95% CL. We also see 
that, in the case of X 1 X 3 production (lower left 
panel) (which is very similar to the cases of X 1 X 4 ) 
X 2 X 3 und x?xj production that we do not show) 


the minimum of the global x^ function for the 
threshold lies between 400 GeV and 1000 GeV, 
again with no upper limit at the 95% CL. It 
should be noted, however, that the optional anti- 
tachyon cut would indeed yield upper limits at 
the 95% CL for those production modes. Refer¬ 
ring back to the bottom right panel of Fig. [T^ 
and the right panel of Fig. we see that slep¬ 
ton pair-production thresholds may well also lie 
below 1000 GeV. In all cases, the expected loca¬ 
tions of the thresholds in the CMSSM, NUHMl 
and NUHM2 are at much higher centre-of-mass 
energies. 

Thus, the accessibility of supersymmetric par¬ 
ticles at e+e“ colliders is vastly different in the 
pMSSMlO and similar non-GUT models, as com¬ 
pared to the simplest GUT-based models. 


7. Conclusions 


We have performed in this paper the first global 
likelihood analysis of the pMSSM using a frequen- 
tist approach that includes comprehensive treat¬ 
ments of the LHC 8 constraints. This analysis re¬ 
quired many developments and extensions of the 
MasterCode framework that are described in ear¬ 
lier sections of the paper. For example, in or¬ 
der to interpret the searches for coloured sparti- 
cles via jets -I- A -I- If^T signatures at LHC 8 , we 
combine searches sensitive to a variety of differ¬ 
ent cascade channels, whose relative probabilities 
depend on other model parameters. By combin¬ 
ing a sufficiently complete set of channels 18 , we 


capture essentially all the relevant decay chan¬ 
nels, and so achieve a reliable LHC 8 coi con¬ 
straint. In the cases of the LHC8ewk constraints 
from searches for electroweak gauginos, Higgsi- 
nos and leptons, we constructed computationally- 
efhcient models for their contributions to the 
global likelihood function that mimic closely the 
more computationally-intensive results from the 
Atom code. A similar procedure was used for the 
LHC8stop constraints from searches for models 
with compressed stop spectra, with the addition 
that we constructed the likelihoods for some sim¬ 
plified model searches using the Scorpion code. 
These procedures have all been validated exten¬ 
sively, as described in the text. 
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Figure 29. The one-dimensional profile likelihood functions for various thresholds in e'^e~ annihilation. 
Upper left panel: The threshold for XiXi production. Upper right panel: The threshold for associated 
X 1 X 2 Production. Lower left panel: The threshold for associated XiX^ production. Lower right panel: The 
threshold for x^ff[ production. 


The results of our analysis of the pMSSMlO are 
described in Section 3, where we provide many 
details of the global likelihood function. We give 
there the parameters of our best-fit pMSSMlO 
point, while cautioning that its squark and gluino 
mass parameters are poorly constrained. On the 
other hand, some of the pMSSMlO parameters in 
the electroweak sector are relatively tightly con¬ 
strained. For example, we find relatively narrow 
ranges of Xj* and slepton masses, which are quite 
light, and that m^o ~ m^o ~ rn.^± in the region of 


parameter space that is preferred at the 68% CL. 
The light spectrum of electroweakly-interacting 
sparticles is preferred by the {g — 2)^ constraint, 
and the neutralino and chargino mass degenera¬ 
cies are then required to obtain a satisfactory cold 
dark matter density. In addition to the best-fit 
point, we have presented and analyzed several al¬ 
ternative pMSSMlO points with low stop, squark 
and gluino masses that may serve as benchmarks 
for LHC Run 2 analyses 

^'*SLHA files [ 37 ] for these points can be downloaded from 
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One of the most striking features of our anal¬ 
ysis is that the pMSSMlO can provide an ex¬ 
cellent fit to {g — 2 )^ while respecting all the 
LHC 8 constraints, something that is not possi¬ 
ble in models with universal soft supersymmetry¬ 
breaking terms at the GUT scale, such as the 
CMSSM, NUHMl and NUHM2. A corollary is 
that there are interesting prospects for exploring 
the preferred region of the pMSSMlO parameter 
space in future experiments. For example, LHC 
searches at 14 TeV have excellent prospects for 
exploring the preferred regions of rriq and rrig, as 
well as light U, e and fl masses. Looking fur¬ 
ther ahead, the {g — 2 )^-friendly regions of the 
pMSSMlO could be explored in detail with an 
e+e” collider operating at 500 to 1000 GeV in 
the centre of mass. In particular, such a machine 
would have a significant discovery potential in the 
preferred region for the lightest neutralino and 
chargino, while those states would be difficult to 
access at the LHC with the searches discussed in 
this paper. Also, we recall that the region of the 
pMSSMlO parameter space that is favoured at 
the 68 % CL after implementing the LHC 8 con¬ 
straints yields relatively large values of cr®^ that 
should be accessible to forthcoming experiments: 
see the right panel of Fig. 

It is a characteristic of the pMSSM that the 
possibility of extrapolation to high renormalisa¬ 
tion scales is not enforced, and indeed we find that 
most of our pMSSMlO parameter sets yield some 
tachyonic sfermion masses at high renormalisa¬ 
tion scales. It is not clear that such models should 
be rejected out of hand [^, but it is reassuring 
that many features of our pMSSMlO fit would, 
nevertheless, be preserved if one required the ab¬ 
sence of tachyons. On the other hand, the pre¬ 
ferred region of the pMSSMlO parameter space 
has non-universal gaugino and sfermion masses. 
The former arise from the tension between (g—2)^ 
(which favours small Mi^ 2 ) and the LHCScoi con¬ 
straint (which favours larger M 3 ) as well as the 
dark matter constraint (which favours Mi ~ M 2 
at the electroweak scale, not at the GUT scale). 
In parallel, sfermion mass non-universality also 
arises from the tension between {g — 2 )^ (which 


the MasterCode website [25]. 


favours small mfi) and the LHC8coi constraint 
(which favours large squark masses). 

It would be desirable to extend our approach to 
more general variants of the pMSSM with fewer 
restrictions on the parameters. For example, it 
would be interesting to relax the assumption of a 
single slepton mass scale: this is unlikely to alter 
the preferred range of the but would have 

important repercussions for dark matter density 
calculations. It would also be desirable to re¬ 
visit in more general pMSSM scenarios the prefer¬ 
ences we have found for neutralino and chargino 
mass degeneracies, and the constraints we find 
in the (M^,tan/3) plane, which are largely in¬ 
direct (being due to the interplay between con¬ 
straints whose combination may have different 
implications in more general pMSSM scenarios). 
However, we think that many features of our 
pMSSMlO analysis would persist in more general 
scenarios. 

Finally, when interpreting the impacts of ex¬ 
perimental searches in our preferred pMSSMlO 
region, it is important to take into account decay 
chains involving an intermediate chargino, which 
is required to be light in order to fulfil the relic 
density constraint. In a large part of our preferred 
parameter space the chargino is almost mass de¬ 
generate with Xij E^nd there are also regions with 
a sizeable mass difference that exhibit distinctive 
decay chains. Therefore, the pMSSMlO moti¬ 
vates interpreting searches not only in terms of 
the minimal decay chains of the simplified mod¬ 
els presently being considered, but also with the 
xf (and possibly also the x®) incorporated in the 
spectrum over a range of low masses. 

We await with interest the verdict of future 
runs of the LHC. 
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